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Abstract 


Strongly  nonlinear  dispersive  waves  described  by  a  general  Klein-Gordon 
equation  with  slowly  varying  coefficients  and  a  dissipative  perturbation  are 
considered.  The  method  of  multiple  scales  shows  that  the  equation  for  wave  action 
and  the  dispersion  relation  (based  on  the  averaged  energy)  are  valid  to  a  higher 
order  than  anticipated.  Thus,  higher  order  terms  are  calculated  from  the  first 
variation  of  wave  action.  The  spatial  and  temporal  slow  modulations  of  the 
phase  shift  are  shown  to  be  governed  by  a  new  equation,  which  is  universal  for 
small  but  otherwise  arbitrary  dissipation.  This  result  extends  to  nonlinear 
partial  differential  equations  the  quite  recent  work  by  the  authors  on  new 
adiabatic  invariants  for  nonlinear  oscillations  governed  by  ordinary  differential 
equations. 
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1.  Introduction, 


Oscillatory  solutions  of  strongly  nonlinear  dispersive  waves  are  quite 
common,  perhaps  the  most  well-known  being  the  cnoidal  waves  for  the  Korteweg- 
deVries  equation.  In  order  to  understand  the  effects  of  a  slowly  varying  medium,  Luke 
[1 ]  in  1966  utilized  the  method  of  multiple  scales  to  analyze  a  model  nonlinear  problem, 
the  nonlinear  Klein-Gordon  equation.  This  extended  Kuzmak's  work  [2]  on  oscillatory 
solutions  of  nonlinear  ordinary  differential  equations  to  the  case  of  nonlinear 
partial  differential  equations.  In  this  way,  the  amplitude-dependent  dispersion 
relation  could  be  used  to  show  that  the  slow  time  and  spatial  evolution  of  the 
amplitude  of  strongly  nonlinear  dispersive  waves  was  governed  by  conservation 
of  wave  action,  the  appropriate  generalization  of  the  adiabatic  invariant  of 
nonlinear  ordinary  differential  equations.  With  a  small  perturbation  representing 
damping,  Whitham  [3]  showed  how  wave  action  was  dissipated. 

Although  the  modulated  phase  shift  for  a  nonlinear  wave  is  part  of  the 
leading  order  long  time  slowly  varying  solution,  it  has  eluded  previous  efforts 
until  now.  For  unperturbed  non-dissipative  systems,  Luke  [1]  correctly  observed 
that  the  phase  shift  satisfied  a  second-order  equation,  one  solution  of  which 
is  a  constant.  However,  the  equation  was  not  stated,  which  apparently  has  led  to 
some  common  misconceptions  concerning  the  phase  shift.  The  same  non-dissipative 
situation  was  reexamined  by  Dobrokhotov  and  Maslov  [4].  They  introduced  some 
imaginative  ideas  that  enabled  the  perturbation  method  to  be  carried  out  to  the 
higher  order  necessary  for  the  calculation  of  the  phase  shift.  Unfortunately, 
their  technique  is  not  as  easy  to  implement  as  ours,  which  we  believe  is  the  reason 
they  incorrectly  omit  a  term  in  their  analysis.  For  the  case  of  strongly  nonlinear 
ordinary  differential  equations  with  small  arbitrary  damping,  we  (Bourland  and 
Haberman  [5])  quite  recently  obtained  the  second-order  differential  equation  for 
the  modulated  phase  shift.  We  showed  that  the  phase  shift  4>(T)  is  determined 
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from  a  new  and  unexpected  adiabatic  invariant  of  universal  form,  independent  of 
the  type  of  small  (but  arbitrary)  damping;  the  small  frequency  modulation  d<f>/dT 
was  shown  to  be  a  constant  multiple  of  the  derivative  (with  respect  to  energy) 
of  the  amplitude-dependent  frequency  of  the  nonlinear  oscillator. 

In  §2,  we  apply  the  method  of  multiple  scales  to  the  nonlinear  Klein-Gordon 
equation  with  slowly  varying  coefficients  and  a  small  arbitrary  perturbation. 
There  we  briefly  rederive  the  well-known  equation  for  the  dissipation  of  wave 
action: 


^  (-u.1)  -  c2(X,T)V.  (kJ)  =  -  R  ,  (1.1) 

where  k  is  the  wave  number  vector,  «  the  frequency,  -ujJ  the  wave  action, 

p 

-c  Jk  the  flux  of  wave  action,  and  R  the  dissipation  of  wave  action.  In  §3, 
we  begin  to  analyze  the  perturbation  expansion  to  a  higher  order.  If  the 
leading  order  perturbation  is  dissipative,  then  we  show  that  wave  action 
satisfies  (1.1)  to  at  least  one  higher  order  than  would  have  been  anticipated. 
Thus,  we  show  that  the  equation  needed  to  describe  the  modulations  of  the  phase 
shift  4>(X,T)  may  be  obtained  in  a  simpler  way,  by  just  considering  the  first 
variations  of  the  well-known  equation  for  the  wave  action.  We  conjecture  that 
this  is  a  general  principle. 

We  obtain  the  phase  shift  in  §4  by  eliminating  secular  terms  at  a  higher 
order  in  our  perturbation  expansion.  Here,  we  show  that  the  leading  order 
perturbation  u^  of  the  leading  order  solution  Ug  satisfies 

Ut  =  Duq  ,  (1.2) 

even 


where  in  §4  we  obtain  D  and  show  it  to  be  a  linear  first-order  differential 
operator  in  the  slow  parameters,  energy  and  wave  number  vector.  We  choose 

a 

D  in  the  appropriate  way  so  that  it  is  easy  to  derive  that  the  partial  differential 

equation  for  the  modulated  phase  shift  <j>(X,T)  is 

*** 

D-j- ( -axJ )  -  c2(X,T)Dy  •  (kj)  =  0  .  (1.3) 

A 

Later  (in  §6)  we  show  that  D  is  the  Taylor  series  operator  in  the  parameters 
E  and  k,  so  that 

° =  Ei  lir  *  V  7k  •  <’-4> 

where  k-j  =  Vcp  is  the  perturbation  of  the  wave  number  and  E-j  is  the  average  of 

A  A 

the  perturbation  of  the  energy.  In  (1.3)  D-j.  and  are  partial  derivatives  of 

A  A  A 

the  operator  D.  (In  one-spatial  dimension  =  D^,  while  in  higher  spatial 

A  A 

dimensions  is  defined  in  §4.)  It  is  interesting  to  note  that  once  D  is 
determined,  (1.3)  can  be  obtained  from  the  wave  action  equation  (1,1)  by 

r\  A  A 

replacing  by  Dj,  V  •  by  Dv«,  and  omitting  the  dissipation.  Equation  (1.3) 

A 

is  a  partial  differential  equation  since  the  coefficients  for  D  depend  on 
4>T  and  V<j>. 

In  §5  we  obtain  equivalent  results  for  the  Klein-Gordon  equation  by 
elementary  (but  innovative)  considerations  of  the  first  variation  of  wave  action 
(1.1)  with  respect  to  perturbations  of  the  wave  number  and  averaged  energy.  We 
show  that  the  perturbation  of  the  frequency  is  accurately  determined  from  a 


higher-order  dispersion  relation  which  includes  the  averaged  perturbed  energy. 
From  these  ideas,  we  obtain 


E1  ^  +  v<t>T  *  VR(- 03  J )  -  c2VE]  •  ^  (kJ)  -  c2JV2c()  =  0  , 

0.5) 

where  E-j  =  (-<£_-  is  the  average  of  the  perturbation  of  the  energy. 

Equation  (1.5)  is  shown  to  be  equivalent  to  (1.3).  Its  form  is  universal ,  since 
it  is  independent  of  the  dissipation.  We  conjecture  that  (1.5)  is  a  new  general 
principle  for  the  determination  of  the  modulations  of  the  phase  shift.  In  the 
case  of  ordinary  differential  equations,  our  [5]  recently  proposed  adiabatic 
invariant  follows  directly  from  (1.5),  and  §5  dramatically  simplifies  its 
interpretation.  Stability  of  these  nonlinear  periodic  waves  are  discussed  in  §7. 

2.  Nonlinear  Dispersive  Waves. 

In  this  section,  we  will  derive  in  the  standard  way  the  equation  for  wave 
action  for  the  perturbed  nonlinear  Klein-Gordon  equation  with  slowly  varying 
coefficients 

utt  "  +  Vu^u,i$’T)  +  e  h(u,ut,Vu,X,T)  +  e2h1  (u,ut,Vu,X,T)  =  0 

(2.1) 

where  V(u,X,T)  is  a  nonlinear  Dotential  depending  on  the  slow-time  scale 

2  N  a2 

T  =  ct  and  slow  spatial  scales  X  =  ex.  The  notation  =  l  -2-*-  is  being  used. 

i=l  3xf 

We  assume  the  potential  admits  oscillatory  solutions  of  (2.1)  (for  all  fixed  X,  T) 
on  the  fast  time  scale  t  without  the  perturbation  (i.e.,  when  h  =  h^  =  0).  Unlike 
Luke  [1]  and  Dobrokhotov  and  Maslov  [4],  we  allow  an  arbitrary  slowly  varying 
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small  0(e)  perturbation,  eh.  To  show  the  effects  of  an  even  smaller  perturbation, 

2 

we  include  the  e  h^  term.  For  the  early  parts  of  this  paper,  we  make  no  further 
assumptions  concerning  these  perturbations.  However,  the  most  interesting  conclusions 
(that  we  will  reach  in  later  sections)  occur  when  eh  is  a  dissipative  perturbation, 
so  that  h  is  odd  in  the  combined  arguments  ut  and  Vu 

p 

[i.e.,  h(u, -ut, -Vu,  X,  T)  =  -  h(u,  u^,,  Vu,  X,  T)]  and  e  h^  is  non-dissipative 
(even  in  the  combined  arguments). 

We  use  the  method  of  multiple  scales  [6]  with  the  fast  scale  and  the  slow 
scales  T  and  X: 

6(X,T) 

¥  =  +  O(X.T) 

T  =  et 

X  =  £  X  , 

and  follow  the  procedure  we  [5]  used  for  ordinary  differential  equations.  This 
form  of  the  method  of  multiple  scales  is  particularly  suited  for  the  eventual 
first  determination  of  the  phase  shift  $(X,T).  For  example, 

ut  =  (0T  +  e*T)lV  +  eUT  ’ 
so  that  (2.1)  becomes  exactly 

(eT  +  £4>j)2uw  +  €[(eTT  +  ^TT)u^  +  2(eT  +  e^Ju^] 

-  c2(V0  +  eV<p)2u^  -  ec2[(V20  +  eV2^  +  2(V0  +  eV<f>)  •  Vu^] 

+  e2uTT  -  e2c2V2u  +  V  (u,X,T)  +  eh  +  e2h]  =  0  ,  (2.2) 


-  6  - 


i « 

where  from  here  on  V  e  J  refers  to  differentiation  of  a  quantity 

with  respect  to  the  slow  spatial  variables.  We  now  introduce  the  perturbation 


expansion 


u(iNX,T)  =  u0(»M,T)  +  eu^.X.T)  +  e2  u^.X.T)  +  ... 


It  is  convenient  to  note  that 


vu(u,X,T)  =  Vu(u0,X,T)  +  (eU]  +  e2u2)  Vuu(u0,X,T) 
2 


+  £  "  Vuuu(u0’*‘T>  +  °<e  ) 


h(u,ut,Vxu,X,T)  =  h(u,eTu^  +  £4>tu^  +  euT  ,  V6u^  +  eV<})U  +  eVu  ,  X,T) 


h ( Uq * ©tUq  » V0u ,X,T) 


+  e[u1hu  +  (eTu,  +  d>Tu0  +  u  )h  +  (veu,  +  V0u  +  Vun)*h  ] 

<p  ip  T  V  ip  u  ~v 


+  o{e)  , 


where  hu  e  hu(uQ,eTu0  ,  uQ  V8,X,T)  ,  hy  e  hy  (u  ,6  u  ,u0  V6,X,T)  , 

ip  ip  t  ip  ip 


Uv  "  X  3(3u/3x.)  @i  h^UO,0TUO. ,U0  V9»i’T^ 

»  *  1  I  Ip  Jp 


Substituting  the  perturbation  expansion  into  (2.2)  yields,  to  leading  order, 
the  nonlinear  oscillator  equation: 
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[6?  -  c2(X,T)(ve)?]un  +  V  (un,X,T)  =  0 


(2.3) 


The  higher  order  terms  become 


L(u1 )=R .  for  i  =  1,2,  ... 


(2.4) 


where  L  =  [eT  -  c2(X,T) (V8)2]  — ^  +  ^uu^O**’1^  1S  t*ie  linearization  of  the 
nonlinear  operator  in  (2.3).  It  is  not  difficult  to  obtain  the  right-hand 
sides.  For  example. 


c2(X,T)V2G]u, 


(2.5) 


-  2[cf>TeT  -  c(x,T)ve*v4>]u0 

w 

-  2[0t  -It  -  c2(X,T)V0-V]uo  -  [0TT 

ip 

-  h(uo,0TuQ  , uq  V0,X,T) 

<1*  ip 


R2  =  -  2[<Pt0t  -  c2(X,T)V0-V<()]u1  -  2[0t  -|j,  -  c2(X,T)V0-V]u. 

-  [0TT  -  c2(X,T)V20]u1  -  [(<j)T)2  -  c2(X,T)(V<J>)2]u0 

^  ~  w 

-  2[<J>t  -  c2(X,T)V4>-V]u0  -  [<j»TT  -  c2(X,T)V2<J)]u0  (2.6) 

ip  ip 

-  ]  Vuuu<V*>T>ul  -  [  ~Z  ■  c2(X,T)V2]„0 

o  I 

-  h1  -  u]hu  -  [0Tu.  +  0Tu  +  uQ  ]h  -  [VGu-j  +  V0uQ  +  VuQ]  •  hy  , 

ip  T  ^ 

where  h-j  is  evaluated  at  the  leading  order. 

The  leading  order  equation  (2.3)  is  a  nonlinear  ordinary  differential  equation 
in  the  fast  variable  \p.  Multiplying  by  uQ  and  integrating  once  yields  the 
conservation  of  energy  equation, 


(2.7) 


\  [6j  -  c2(V0)2](uo  )2  +  V(uq,X,T)  =  E(X.T)  , 


where  E(X,T)  satisfies  a  well-known  equation,  which  we  will  rederive.  Standard 

phase  plane  analysis  shows  that  the  solution  Uq(i/j ,X ,T )  is  periodic  in  <p 

for  appropriate  potentials  V(Uq,X,T),  oscillating  between  Uq  an(j  Uq  with 

mi  n  max 

V(uQ  »X,T)  =  E  and  a  similar  expression  involving  uQ  .  As  Luke  [1]  showed 
umi n  ~  umax 

and  we  [5]  repeated,  u n(i^,X,T)  can  be  defined  as  an  even  function  of  ip  with 

~  3u0 

Uq(0,X,T)  =  uQ  or  Uq  and  -r—  (0,X,T)  =  0.  The  period  of  the  oscillation 
mi  n  max  dv 

in  the  fast  variable  is  constant  [1],  and  we  normalize  it  to  1  (although  2tt  is 


often  used)  which  gives 


(eT)2  -  c2(x,T)(ve)2  =  2 

1  -  J, 


un  (x,t) 


[Uq  (x,t)  /2 
min  ~ 


~  [E(X,T)  -  V(u0,X,T)]"2du0 


(2.8) 


I 


to 


In  this  situation,  the  notation  k  =  V9  for  the  wave  number  vector  and  w  =  -Q-j- 
for  the  frequency  is  appropriate.*  Equation  (2.8)  may  be 

regarded  as  the  dispersion  relation,  w  =  u(k(X,T) ,E(X,T) ,X,T) ,  for  the  propagation 

of  nonlinear  waves  in  a  nonuniform  media  as  described  by  (2.1). 

Some  well-known  results  on  the  linearized  operator  are  needed  to  analyze 

the  higher  order  equations,  L(y1-)  =  R. .  Taking  the  derivative  of  (2.3)  with  respect 

to  ip  yields  L(uq  )  =  0,  showing  that  Uq  is  a  homogeneous  solution.  Since 

un(<P,X,T)  is  an  even  periodic  function,  un  is  an  odd  periodic  homogeneous  solution; 
0  ~ 

usually  it  is  the  only  periodic  solution  [1],  In  order  for  the  solutions  of  (2.4) 
to  be  periodic,  the  Fredholm  alternative  implies  that 


In  comparing  this  to  the  case  of  ordinary  differential  equations,  note  that  aJ0<jie=“u>pde* 


I 


v 

I# 

I 

£ 
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fl 


uA  R,  dv  =  0 


n  0,1 

0  ip 


(2.9) 


equivalent  to  the  usual  secularity  condition. 
For  i  =  1 ,  using  (2.5),  (2.9)  becomes 


1 


-20TeT  -  c  (X,T)V6'V<J>]  u0  uQ  dip 

■T)  ipip  ip 

-2[0J  fy  -  c2(X,T)ve*V]  uQ  UQ  dip 

j0  ip  ip 

-[eTT  -  c2(x,T)v2e]  j  un  un  dip 
1  1  ~  JO 

rl 


h(Ug, 0-pUg  ,Uq  V@,X,T)dtp  -  0 

0  ip  ip  ~ 


The  first  integral  is  zero.  The  next  two  integrals  combine  to  yield  Whitham's 
[3]  well-known  partial  differential  equation  describing  the  possible  dissipation 
of  wave  action: 


jr  [eT  (un  )2d<p]  -  c2(x,t)v.  [(ve)  f\un  )2dip] 
jn  u,i,  -  j  o  ^ 


0  > 
,1 
'0 


(2.10) 


h(u0,eru0  »  Ug  V0,X,T)Uq  dip  =  0  . 

^  V  Ip 


We  define  the  latter  integral  to  be  the  dissipation  R  and 


J  = 


rl 


(un  )  dip 

0  uip 


(2.11) 


so  that  (2.10)  may  be  written  as 


1 


fjr  (eTJ)  -C  V  •  (JV6)  +  h(u0,eTu0  ,u0  ve.x.T)  dip  =  o  .  (2.12) 

•*0  ip  ip 


9 


a 


Si 


a 


>2 


9 


1 


% 


r *  v 
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For  ordinary  differential  equations  [2],  this  simply  reduces  to  the  possible 
dissipation  of  action  (an  adiabatic  invariant  if  h  =  0). 

3.  Higher  Order  Wave  Action. 

In  this  section,  we  show  that  wave  action  satisfies  ( 2 . 1 0 ) to  at  least 

one  higher  order.  We  will  use  this  in  §4  to  determine  an  equation  for  the 

modulations  of  the  phase  shift  <J>( X ,T ) .  To  obtain  these  results,  we  eliminate 

** 

2 

the  secular  terms  in  the  0(e  )  equation: 
fl 

1  Rr,  Un  dtp  =  0  ,  (3.1) 

J0  Uip 


3 


where  is  given  in  (2.6).  As  with  the  case  of  ordinary  differential 

equations  ([4], [5]),  the  following  difficult  term  can  be  calculated  by  integration 

by  parts,  then  by  recalling  the  definition  of  L,  and  finally  by  using  L(u-j)  =  R1 : 


1  ?  f1 

.Vuuu(V~X>T)ul  “o*  ■  -2  (.W*>T)ul“l1d* 

)  ip  J0  \p 

_  _  2  ^  n  ' ••  '  ' .  .2  _2,.2 


CL ( u -j )  -  (go  -c  k  )U|  ]u,  d^ 


=  2  [  (2(d)TeT  -  c2ve-V0)uo  +  2(eT|  -  c2ve*v)uQ 

■'0  ip\p  *  [p 

2  2 
~  ^  n4-  r 


+  (Qyy  "  c  V  0)Uq  +  h(uQj0jUQ  »  Uq  V0,X,T)}  dip 

$  >  ^  ~  \p 


where  the  periodicity  of  some  expressions  has  also  been  used.  In  this  manner, 
from  (.3.1),  we  obtain 


W’J  "J  Ml  ■_»  rjrrcr  s  -j-  v-  rw. 


■§f  [4>t  (uq  f  dip]  -  c2V*[V4>  (uq  ) 2  d  iJj] 

J0  \p  JO  ip 

2  JT  [eT  |'U0IU11<)*]  •  c2v-t76  |'u0,ul,d,)’] 

J0  ip  \p  1,1  i|j 


'0  >  > 


f1  r1  f1 

+  nhuld*  +  nulVhud*  +  L(6Tu1,  +  Vo,  +  uo>  uo(hv 

J0  ip  J0  \l>  JU  U)  [I;  T  ill 


(3.2) 


+  (veu.  +  V<}>u0  +VuJ  uq  *h„  dijj  +  h^Q  d^  =  0 

JQ  ip  i|>  ip  JQ  \p 


A  similar  calculation  was  presented  by  Dobrokhotov  and  Maslov  [4],  They 
incorrectly  calculated  one  term  and  thus  did  not  obtain  (3.2).  Furthermore, 
their  results  were  for  unperturbed  systems  (h  =  h^  =0).  Equation  (3.2)  is 
not  a  difficult  equation;  by  comparing  it  to  (2.10),  we  see  that  most  of  the 
terms  are  the  first  variation  of  the  wave  action  equation  (2.10),  as  we  now 
show.  Wave  action  is  a  concept  valid  for  the  leading  order  terms  in  a 
perturbation  expansion.  The  wave  action  defined  in  terms  of  total  frequency, 
total  wave  number  vector,  and  total  u  will  not  satisfy  (2.10).  However,  we 
expect  there  might  be  an  0(c)  correction: 


Yf  [(0T+  e<i>T)  K  +  eu1  +  . . . )  d 4*] 
•'0 


-  c2v*[(ve  +  e7<f>)  f  (un  +  eU.  +  . . .  )2  d^] 
J0  v  > 


+  h(uQ  +  eu^  +  . . . ,  (6j  +  e4)j)  (Uq  +  eu-|  +  ...),  (V6  +  cV<{>)(Uq  +  eu-|  +.  .),X,T)cfy 

Jq  ^  ip  p  ip 


=  eP  . 
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We  compare  the  above  to  (3.2)  using  (2.10)  and  determine  that 


P  =  - 


J0 


(u 


0  \  +  VU0  *  hV  +  h1^u0  + 

T 


This  gives  us  a  simple  method  to  obtain  higher  order  corrections  to  wave  action. 

In  general,  action  is  not  satisfied  to  higher  order  (i.e.,  P  f  0).  However, 

2 

if  eh  is  a  dissipative  perturbation  and  e  h-j  is  not  a  dissipative  perturbation, 
then  the  integrals  above  vanish  (since  hv>  h^  and  h-j  are  even).  In  this  important 
case,  our  perturbation  method  has  proved  that  wave  action  is  dissipated  in  the 
same  way  as  (2.10)  to  at  least  one  higher  order  than  would  have  been  anticipated. 


4.  Modulated  Phase  Shift. 


We  restrict  our  attention  to  the  situation  in  which  eh  is  dissipative  (and 


2 

thus  odd  in  i^)  and  e  h^  is  non-dissipative  (and  thus  even  in  4>).  To  obtain 


an  equation  for  the  modulations  of  the  phase  shift  <J>(X,T),  we  determine  u-j  and 


substitute  it  into  (3.2).  We  will  obtain  a  universal  form.  Because  Ug  is  even 


in  \i>  and  h  is  odd  in  \p,  we  only  need  the  even  part  of  u,,  u.  .  From  (2.4)  and  (2.5) 

even 


L(u 


)  =  -  2(-  uxf>T  -  c  k  •  V<f>)ur 


(4.1) 


even 


W 


We  can  easily  solve  (4.1)  by  noting  the  result  of  differentiating  (2.3)  with 
respect  to  E  (keeping  X,  T,  and  k  fixed): 


L(u  )  =  -2WWe  u 
E  ^ 


(4.2) 
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« 


since  w(E,k,X,T).  This  idea  was  used  by  Luke  [1]  and  Dobrokhotov  and  Maslov  [4]. 

However,  we  can  also  take  partial  derivatives  with  respect  to  each  component  of  the 

2 

wave  number  vector  k.  Since  wV.co  =  c  k  from  (2.8),  it  follows  that 

K  *v 


L(Vku0)  =  -  2  (u)Vkoj  -  c  k)uQ  =  0 


(4.3) 


Thus,  Uq  is  independent  of  k  (as  is  also  clear  from  (2.7)  and  (2.8))  since 
VkuQ  is  even  and  periodic.  [The  only  even  homogeneous  solution  is  not  periodic 
(see  [1]  or  [5]).]  By  comparing  (4.1)  with  (4.2): 


v  • 7, 


Thus,  u1  can  be  obtained  by  operating  on  uQ: 
even 


u,  =  Du0 
even 


(4.4) 


(4.5) 


It  is  convenient  and  allowable  (but  not  necessary)  to  include  derivatives  with 

A 

respect  to  k  since  VkUQ  =  0.  Therefore,  from  (4.4),  we  define  D  as  the 
following  linear  first-order  differential  operator  in  the  slow  parameters  E  and  k: 


A  "^T  ”  Vkw  *  ^  3 


#+rvk  * 


(4.6) 


where  y  is  arbitrary  since  VkUg  =  0.  However,  we  will  discover  that  our 
derivation  is  simplified  for  a  specific  choice  in  which  y  f  0. 


We  note  that  for  any  linear  first-order  differential  operator 


h(u0,-uiu0  ,ku_  ,X,T)uq  dip  = 
rl  rl 

hu-|  dip  +  uQ  [h  u,  +  h  (-uu,  -  D(w)u0)+ (ku1  +  0(k)uQ  )  •  hy3dip 

;0  ^  ;0  ij;  ip  \p  ~  V  ~  V  ~ 

(4.7) 


since  derivatives  in  the  fast  variable  ^  commute  with  D  and  since  only  the 
even  part  of  u-j  contributes  to  the  integrals.  These  are  precisely  the  terms 
that  appear  in  (3.2)  if 


D(uj)  =  -<pj 

B(k)  =  v*  , 


(4.8a) 

(4.8b) 


which  will  greatly  simplify  our  derivation.  Thus,  we  will  show  only  one  y 
can  be  chosen  to  satisfy  (4.8).  From  (4.6) 

A 

D(u)  =  -<tj  -  Vko)  •  V<p  +  y  *  V^uj 
D(k)  =  y  . 

Thus,  (4.8b)  is  satisfied  only  if 

y  =  V<^  , 


which  fortunately  does  Sc'iSfy  (4.8a),  so  that 
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-  -<pj  -  V.  a)  •  V<j) 

0  =  ~ -q -  •  <«.9> 

(In  §6,  we  show  that  D  is  the  Taylor  series  operator  in  the  parameters  E  and  k. 

A 


) 


Substituting  u-j  = 
*  even 
because  D  is  a  linear 


DUq  into  (3.2)  yields  relatively  easy  calculations 
first-order  differential  operator.  For  example. 


A  f  1  2  1 

D  (un  )  d^  =  2  un  u.  dijj  . 
;0  u\p  J0  Uip  V 


Thus,  (3.2)  simplifies: 


jf  (4>tJ)  -  c2V  •  (V4kJ) 

+  (0-j-DJ)  -  c2V  * 

where  again  J  is  given  by  (2.11). 


A  [ 

(VQDJ)  +  D  hun  dip 

J0  Uip 

However,  from  (4.8) 


0  , 


D(0tO)  =  0-j-DO  +  D(0t)J  =  0tDJ  +  <pjj 

A  A  A  A 

D(V0J)  =  V9DJ  +  D(ve)j  =  V8DJ  +  VcfrO  , 


(4.10) 


since  again  D  is  a  first-order  differential  operator.  Thus,  the  conservation 
form  of  the  equation  becomes 


D(0tJ)  -  c2V-D(V0J)  +  D 


1 

hun  dip 

J0  > 


0 


(4.11) 


This  is  the  partial  differential  equation  for  the  modulated  phase  q>( X ,T)  since 

A 

D  is  given  by  (4.9).  This  result  is  particularly  pleasing; it  can  be  obtained 


from  the  well-known  dissipation  of  wave  action,  (2.10)  or  (2.12),  by  replacing 
~  by  0,  V  •  by  V  •  0  and  the  dissipation  integral  by  D  operating  on  it. 

A 

However,  (5.11)  can  be  further  simplified  by  just  operating  D  on  (2.12): 

D  f  hu0  dii>  =  -  D  (eTJ)  +  c2DV*[V6J]  . 

Jo  3T  T 

The  above  can  be  used  to  "eliminate"  the  dissipation  term  in  (4.11): 

^D(6tJ)  -  O^(0TJ)  -  c2V-D(V6J)  +  C2DV-(V9J)  =  0  .  (4.12a) 

Therefore,  we  obtain  a  particularly  simple  representation: 

Dt(6tJ)  -  c2Dv*  (V6J)  =  0  ,  (4.12b) 

A  A  A  A 

since  DT  =  -gy  D  -  D^y  and  Dv  =  V  •  0  -  DV.  Equation  (4.12)  has  been  derived 
for  rather  arbitrary  dissipation;  its  form  is  universal.  The  dependence  on 
dissipation  is  due  to  the  dissipation  of  wave  action  (2.12).  This  result  suggests 
the  following  algorithm  for  obtaining  (4.12):  Starting  from  dissipation  of  wave 

rs  A  A 

action  (2.12),  replace  yy  by  Dj,  V  by  Dy,  and  ignore  the  dissipation  integral. 

/V  A  /A 

Note  that  in  one  spatial  dimension,  Dy  =  Dx,  which  is  as  simple  as  Dy.  Similar 
results  are  valid  for  Korteweg-deVries  type  equations  [7].  For  clarity,  we  note 
that  after  using  (4.9)  and  V^J  =  0  [see  (7.6)]  the  equation  for  the  modulated 
phase  shift  4>(X,T)  becomes 


J 
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-4>T  -  VkU3  •  V<j) 

0  =  ( - - - )j  (-  co  J)  +  V<J>-|-  •  V^(-  coJ) 


OOr- 

c 


?  -4>t  ■  •  v<j)  . 

-  Cz  v  - )  •  £  (kJ) 


Wr 


2  2 

-  c  J  V  <j>  . 


(4.13) 


Here,  for  example,  -^r  refers  to  keeping  only  X  fixed  (not  E  or  k). 

2 

In  the  spatially  independent  case  (c  =  0),  (2.1)  is  a  nonlinear  ordinary 
differential  equation  and  (4.13)  becomes 

0  • 

the  adiabatic  invariant  for  "arbitrary"  dissipative  systems  that  we  [5]  recently 
discovered. 


5.  A  Higher-order  Accurate  Dispersion  Relation  Depending  on  Averaged  Energy, 

In  §4,  we  derived  (4.13),  the  equation  for  the  modulations  of  the  phase 
shift  <t>(X,T) .  Here,  we  will  show  how  this  result  can  be  obtained  from  two 
intuitive  physical  principles  (and  justify  these  ideas  mathematically).  We 
have  already  shown  one  of  these  physical  principles:  the  wave  action  equation 

(-u)J)T  -  c2V  •  (k  J)  +  R  =  0  (5.1) 

is  valid  to  at  least  one  higher  order  in  an  asymptotic  expansion  in  powers  of  e, 
if  eh  is  a  dissipative  perturbation  and  e  h|  is  a  non-dissipative  one. 


4 
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We  conjecture  (and  later  verify)  a  second  physical  principle,  one  based 
on  the  averaged  energy.  We  have  defined  via  (2.7)  a  leading  order  slowly 
varying  energy  Eg(X,T)  in  terms  of  the  leading  order  frequency  Wg,  wave 
number  vector  kg,  and  solution  u^.  We  will  first  introduce  more  accurate 
expressions  for  the  energy.  There  are  two  different  methods  for  doing  this. 
We  can  define  a  mathematical  expression  similar  to  (2.7): 

Emath  -  {  C(eT  +  e4>T)2  -  c2(V0  +  eV<|>)2]u2  +  V(u,X,T)  ,  (5.2a) 


based  on  the  complete  frequency,  wave  number  vector,  and  solution.  Alternatively, 
we  can  use  the  physical  expression  for  energy  density: 


physical  "  2  (ut)2  "  2  c2(vxu)2  +  V(u’*’T)  » 


which,  using  the  multiple  scale  approach,  can  be  rewritten  exactly  as 

2 

^physical  *  7  [(0T  ♦  c«T)u„  *  -  y  [(76  +  e7*)u„  +  eVu]2  a  V(u.X.T)  . 

(5.2b) 

It  is  clear  that  both  equal  Eg  to  leading  order  [defined  by  (2.7)].  It  is 
straightforward  to  show  that  both  Ep(-,ySjcai  and  ^nath*  as  cleir’ned  above» 
oscillate  on  the  fast  scale  with  a  small  0(c)  amplitude.  Moreover,  they  differ 
at  0(e): 


'physical  ”  math 


*  eu. 


(6tuo. 


2 

-  c  ve 


vu0) 


0(e2) 


Both  oscillate  on  the  fast  scale,  so  that  we  might  expect  their  averages 

fl  fl 

to  be  important.  Since  uQ  uQ  dtp  =  0  and  uQ  VuQd^  =  0  by  parity 

■*0  ip  UT  ■’0 


considerations,  we  note  that 


<E  .  .  “  <  E  ,,>  =  0(£  )  * 

physical  math  v  * 


where  <  >  =  dii  represents  averaging  over  the  fast  nonlinear  oscillation. 

Thus,  we  can  use  either  expression  for  the  averaged  energy  and  obtain  identical 
results  valid  to  at  least  one  higher  order  in  powers  of  e. 

By  expanding  either  energy  (5.2a)  or  (5.2b)  to  0(e)  and  by  integrating, 
we  obtain  the  averaged  perturbation  E^: 

Ei  =  {-Vt  '  A)*V<J,)  [ uo2  +  (u)o2  "A? }  U.M*  +  fnVld^ 


J0  >  'Uj 


■'0  0 


(5.3) 


This  expression  can  be  simplified  using  our  results  from  §4.  Only  the  even 

A 

part  of  u,  is  needed  in  (5.3),  namely  from  (4.5)  u,  =  Dufl.  If  we  operate 
/v  even  u 

D  on  (2.7) ,  we  obtain 


Ug-  C  kg2  )Ug  DUq  +  Ug  [a)gD(u)g)  ~  C2kg*D(kg)]  +  V  DUg  =  D  Eg  , 

ip  \p  \p  0 

which,  after  averaging  (integrating)  and  using  (4.8),  yields  an  expression 


for  E1  from  (5.3) : 


E1  "  D  E0  ‘ 


-<|)j  -  v^Wg  *  V4> 


(5.4) 


I 


We  will  show  that  (5.4)  is  an  immediate  consequence  of  a  simple  physical 
conjecture.  The  local  amplitude-dependent  dispersion  relation  (2.8), 


ojq  =  oo[kg(X,T),  Eq(X,T),  X,  T],  is  well-known  to  be  valid  to  leading  order. 
Intuitively,  we  might  expect  that  the  dispersion  relation  is  valid  to  at 
least  one  higher  order  if  both  the  more  accurate  wave  number  vector 
k0  +  e  k-j  (with  k-j  =  v<j>)  and  the  more  accurate  averaged  energy  Eg  +  e  E-j  are 
used.  We  would  then  conjecture  on  the  basis  of  physical  intuition  that  the 
perturbed  frequency  w,  =  -  <f>—  is  obtained  from  the  elementary  Taylor  series  of 
this  highly-accurate  dispersion  relation: 


to. 


=  E, 


U), 


0E  +  ~1 


Vo 


(5.5) 


Since  this  is  equivalent  to  (5.4),  we  have  proved  the  validity  of  a  second 

physical  principle:  the  dispersion  relation  is  valid  to  at  least  one  higher 

order  if  an  accurate  expression  for  averaged  energy  is  utilized  in  (2.7). 

Frequency  changes  are  only  due  to  the  perturbed  wave  number  vector  and  to  the 

averaged  perturbed  energy.  This  is  valid  only  if  the  perturbation  eh  is 
2 

dissipative  and  e  h-j  non-dissipative  (for  otherwise  we  would  expect  further 
frequency  changes  due  to  these). 

We  use  this  result  to  obtain  the  modulations  of  the  phase  shift  by  analyzing 
the  first  variations  of  the  wave  action.  Since  (5.1)  is  valid  to  at  least 
one  higher  order,  we  will  substitute  k  =  kg  +  e  k^  and  to  =  tog  +eto^  into  (5.1), 
where  k-j  =  V<J>  and  =  -  <f>T  .  Furthermore,  we  must  use  the  more  accurate  integral 
R  and  J  defined  in  terms  of  total  k,  w,  and  u.  However,  it  can  be  shown  from 
§4  that  R  and  J  are  functions  of  k  and  the  average  energy: 


1  _t A  .«*-**,  _< 4  a*a  .«*  .•<  .fk  a'|  _  |  -**  .»>  >■»  .»fc  .»!  »  4^|'|  |  |  |  1^'Lft'Ll  «-A  l-l 
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R  -  Rn  +  e  (E-|  Rq  +  k-j  •  V^Rg)  +  0(e C) 


J  =  JQ  +  e  E-j  +  Oiz^)  , 


since  J  is  independent  of  k  [see  (7.6)].  In  this  way,  (5.1)  becomes 


-  [(wg  +  cgjjMJq  +  e  E^  JQ  )]t 


c  v*nkn+ck^)(Jg+  e  E  i  Jq  )] 


(5.6) 


+  Rn  +  e(Ej  Rq  +  kj  ’V^Rg)  +  )  =  0  • 


The  leading  0(1)  terms  cancel,  and 


many  of  the  0(c)  terms  cancel  since  the  partial  derivatives  with  respect 


to  E  and  k  of  the  leading-order  wave  action  equation  are  also  valid.  Thus, 


the  remaining  0(e)  terms  must  be  in  balance: 


E-|  (-w J)^  +  V<t>T  •  Vk(-uJ) 


(5.7) 


-  c2V  E1  •  k  JE  -  c2  J  V2  <p  =  0  , 


where  the  zeroth-order  subscripts  have  been  dropped.  We  call  (5.7)  the 


eouation  for  the  variation  of  wave  action.  When  the  expression  for  averaged 


perturbed  energy  E-j  is  used,  (5.4),  we  obtain  (4.13),  the  linear  partial 
differential  equation  for  the  modulations  of  the  phase  shift.  Since  (5.7)  is 


the  appropriate  first  variation  of  the  wave  action  (5.1),  it  is  in  essence 


HIT  iT*  rors  -  -  -  -,v  rj'X'j » wt  v*  vr  nwwwwwjww.Byw  w  Uggltpmwv  mwirir.  TVfv^?rffr3^^:m^7y' 

\ 
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the  linearization  of  (5.1).  Thus,  the  characteristic  velocities  of  (5.7) 
will  be  the  same  as  those  for  (5.1),  as  shown  in  §7. 

In  summary,  the  equation  for  the  phase  shift  can  be  obtained  from  the 
correct  higher-order  accurate  versions  of  the  dispersion  relation  and  the 
wave  action  equation. 

A 

6.  A  Linearization  Principle  for  D. 

In  §4,  we  noted  that  calculations  needed  for  higher-order  perturbations 

A  A 

were  simplified  if  we  observed  that  u,  =  Du0,  where  D  is  given  by  (4.9). 

even 

Here,  we  will  show  that  this  operator  follows  from  an  elementary  understanding 
of  the  linearization  of  the  solution.  For  the  Klein-Gordon  equation,  to 
leading  order,  nonlinear  oscillatory  waves  satisfy 


Uq  -  u q ( ip ;  E,  k,  X,  T) 


(6.1) 


where  the  usual  dispersion  relation  (2.8)  is  satisfied 


OJ 


=  u(E,k,X,T)  . 


(6.2) 


If  the  perturbation  is  of  the  restricted  class  (dissipative)  discussed 
throughout,  then  higher-order  terms  do  not  effect  the  phase  modulations 
(other  than  via  wave  action).  In  this  case  and  others,  some  of  the  higher- 
order  terms  may  be  obtained  simply  by  perturbing  the  parameters  E  and  k  in 
(6.1)  and  (6.2).  By  this  Taylor-series  approach. 
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=  El  \  +  ~1*  V*<  U0 
even  E 

(6.3) 

3 

o 

• 

-*7 

+ 

LU 

3 

LU 

II 

(6.4) 

Here  E-|  represents  (as  in  §5)  an  averaged  perturbed  energy.  We  again 
introduce  k-j  =  V<f>  and  ,  in  which  case  (6.3)  and  (6.4)  can  be 

combined  to  yield 

A 

u-|  =  D  Ug  ,  (6.5a) 

even 


where 


D  = 


!srvk 


-  V  0  •  tin 

_  j 


Vif,  •  Vk  » 


(6.5b) 


agreeing  with  (4.9).  As  defined  this  way,  it  automatic-51  ly  follows  that 

A  /V  A 

D(k)  =  ki  =  V*.  D(u.')  =  oj-|  =  -  <t>y  ,  and  D(E)  =  E-j 

The  equation  for  the  modulations  of  the  phase  shift  can  be  obtained 
by  higher-order  perturbations  of  the  wave  action  equation  and  dispersion 
relation.  These  principles  are  easily  extended  to  other  systems,  such  as 
KdV  type  equations  [7],  which  has  two  amplitude  parameters.  Furthermore, 
the  equations  for  the  phast  modulations  of  harder  problems,  such  as  nonlinear 
ones  with  multiple  fast  phases,  should  be  analyzed  in  these  ways  as  well. 

If  perturbations  are  not  necessarily  restricted,  then  the  same  type  of 
linearized  system  will  occur  with  nonhomogeneous  terms  (representing,  for 
example,  non-dissipative  effects). 
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7.  Characteristic  Velocities  and  Stability. 

In  this  section,  we  first  show  that  the  characteristic  velocities  for 
wave  action  are  the  same  as  the  characteristic  velocities  for  the  modulations 
of  the  phase  shift.  Then,  we  calculate  these  characteristic  velocities  for 
nonlinear  waves  periodic  in  one-spatial  dimension.  We  show  that  Klein-Gordon 
waves  are  stable  for  hard  potentials  (i.e.,  >  0),  while  wildly  unstable 

for  soft  potentials  (i.e.,  u^/to  <  0).  These  soft  Klein-Gordon  waves  will  be 
shown  to  satisfy  an  elliptic  partial  differential  equation,  and  hence  are  not 
well-posed  as  a  slowly  varying  wave,  in  the  same  manner  in  which  the  Benjamin- 
Feir  instability  arises  for  water  waves  and  certain  nonlinear  Schrodinger 
equations  (see  Whitham  [3]).  For  this  unstable  case,  nonlinear  waves  with 
nearly  one  wave  number  and  frequency  will  not  persist.  Instead,  solitary  waves 
(if  they  exist)  or  multi-phased  spatially  periodic  waves  (if  they  exist)  might 
develop  (Newell  [8]). 

The  wave  number  vector,  frequency,  and  amplitude  of  a  nonlinear  periodic 
wave  slowly  evolves  in  space  and  time  according  to  the  coupled  system  composed 
of  the  dispersion  relation  (2.3),  conservation  of  waves 

kT  +  Vw  =  0  ,  (7.1a) 

and  the  principle  of  wave  action 

(-  ai  J)T  -  c2V  •  (kJ)  +  R  =  0  .  (7.1b) 

By  considering  the  frequency  w  as  a  function  of  k,  E,  X,  T,  via  the  dispersion 


rVCV  V  VA 
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relation,  the  system  (7.1)  is  expressed  in  terms  of  the  fundamental  unknowns 
E  and  k: 


k.j  +  +  k^  •  Vkw  =  “  3  oj/ 3X ^ 


Et(-ox3)e  +  ky  •  Vk(-  ojJ)  -  c2JV  •  k  -  c2(k  •  VE)Je 


R  -  ^  (-oj  J)  +  c^Vx  •  (kJ)  , 


(7.2a) 


(7.2b) 


where  [only  in  the  rhs  of  (7.2)]  ,  gy  ,  and  refer  to  differentiation 

keeping  k  and  E  fixed.  Equation  (7.2b)  has  been  simplified  somewhat  since  J 
does  not  depend  on  k  [see  (7.6)].  By  comparing  equations  for  the  wave  action 
(7.2)  to  the  phase  shift  (4.13),  it  is  apparent  that  the  characteristic  velocities 
are  identical. 

We  determine  the  characteristic  velocities  for  slowly  varying  Klein-Gordon 
waves  in  one  spatial  dimension.  In  this  case,  (7.2)  becomes 


ky  o.'E E^  ^  iok  k  x  l.o. t. 


(7.3a) 


(-oxOyEy  +  (-ujJ)kky  -  c2Jkx  -  c2kJEEx  =  l.o. t.  ,  (7.3b) 


where  l.o.t.  designates  lower  order  terms,  which  are  not  necessary  to  calculate 
the  characteristic  velocities.  Equations  (7.3)  form  a  system  of  quasi-linear 
first-order  partial  differential  equations.  By  the  usual  method  of  diagonal ization 
(described  and  calculated  by  Whitham  [3]),  we  obtain  the  two  characteristic 
velocities 


mm 


m 


ws 


tiv* 


where  ojcj^  =  c^k  has  been  used  to  derive  this  result  as  well  as 


(u>^  -  c^k*')^  +  =  1 


(7.5) 


The  latter  equation,  (7.5),  is  derivable  from  conservation  of  energy. 


A  (u2  -  c2kZ)J  -  f  (E-  V)dii  =  lJ-  c2k2)'5  f  0,nax  JT  (E-  V)*5  du„  , 

2  J°  J“0  . 

min 

(7.6) 

2  2  2 

when  the  dispersion  relation  (2.8)  is  used.  From  (2.8),  ou  -  c  k  only  depends 
on  X,  T  and  E.  Thus,  J  does  not  depend  on  k,  a  result  we  have  found  useful 
occasionally  in  this  paper. 

The  characteristic  velocities,  given  by  (7.4),  determine  the  stability 
and  wel 1-posedness  of  nonlinear  periodic  waves  for  the  Klein-Gordon  equation. 
The  sign  of  is  important,  since  J  >  0.  Note  that  is  determined  from 
the  potential  [see  the  dispersion  relation  (2.8)].  If  uv/w  >  0,  which  we 
define  as  corresponding  to  hard  potentials  since  the  frequency  increases  with 
the  wave  amplitude  or  energy,  then  the  nonlinear  periodic  waves  are  stable. 

The  initial  value  problem  is  well-posed  (see  [5]).  However,  if  w^/u)  <  0, 
corresponding  to  so-called  soft  potentials,  the  velocities  have  a  non-zero 
imaginary  part,  implying  that  the  partial  differential  equations  are  elliptic. 
In  this  case  the  partial  differential  equation  is  not  well-posed  as  an  initial 
value  problem;  the  assumption  of  a  slowly  varying  periodic  wave  is  not  valid. 


The  periodic  wave  would  immediately  be  wildly  unstable,  yielding  some  other 
type  of  solution.  Thus,  the  phase  shift  analysis  in  this  paper  is  applicable 
only  for  hard  potentials  (u^/w  >  0). 
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Abstract 

The  phase  shift  and  corresponding  small  frequency  modulation  for  weakly 
dissipated  nonlinear  oscillators  with  slowly  varying  coefficients  is  calculated 
for  the  first  time.  This  extends  and  corrects  earlier  work  by  Kuzmak,  Luke, 
and  Dobrokhotov  and  Maslov.  A  new  adiabatic  invariant  is  derived,  the  ratio  of 
the  perturbation  of  the  frequency  to  the  derivative  (with  respect  to  energy)  of 
the  amplitude-dependent  frequency  of  the  nonlinear  oscillator. 


1- 


1.  Introiu  ction. 

Kuzmak  [4]  analyzed  the  l^ng  time  effect  of  weak  damping  on  strongly 
nonlinear  slowly  varying  oscillators.  In  his  important  work  the  leading  order 
solution  was  obtained.  However,  it  was  unjustifiably  assumed. that  the  phase 
shift  of  the  oscillations  was  constant,  and  thus, that  the  frequency  was  not 
perturbed.  Since  that  work,  refinements  have  been  suggested  ([1].[2],[3],[5]), 

Luke  [5]  developed  a  more  systematic  approach  to  the  problem  using  the  method  of 
multiple  scales  that  allowed  the  calculation  of  higher  order  terms.  His  paper, 
nowever,  only  investigated  the  use  of  the  method  on  the  nonlinear  Klein-Gordon 
equation  and  did  not  examine  the  dissipative  perturbations  discussed  here.  He 
correctly  observed  that  the  phase  shift  satisfied  a  second-order  differential 
equation,  one  solution  of  which  is  a  constant.  He  did  not,  however,  state  the 
differential  equation  or  the  form  of  the  second  solution.  This  has  led  to  some 
misconceptions  among  some  authors  about  the  phase  shift.  Dobrokhotov  and  Maslov 
[1]  examined  the  same  (non-dissipative)  equation  as  Luke  [5],  Using  some 
excellent  new  ideas,  they  determined  a  differential  equation  for  the  phase  shift. 
Unfortunately,  they  inadvertantly  omitted  a  term  in  their  dissipationless  analysis 
which  caused  an  error  in  their  result.  In  this  paper,  we  will  derive  the  correct 

second-order  differential  equation  for  the  Dhase  shift  and  do  it  in  a  simpler 
manner,  including  arbitrary  damping.  We  also  obtain  the  solutions  to  this  equation 

and  show  how  it  leads  to  a  new  and  unexpected  adiabatic  invariant  for  dissipative  equations. 

In  §2,  we  analyze  second-order  nonlinear  differential  equations  of  the  form 


with  coefficients  varying  on  the  slow  time  T  =  c t.  Here  V(y,T)  is  a  nonlinear 


v'-V-V-  Vv'/v'v'‘V-Vvv’' VkvS-^''^SrS-%N^ 


potential  that  admits  periodic  solutions  when  c  =  0.  The  term  h(y,^-,T)  is 
assumed  to  be  odd  in  dy/dt  in  order  to  represent  an  arbitrary  small  nonlinear 
damping.  This  generalizes  the  approach  of  Kuzmak  [4]  who  only  included  quasi- 
linear  damping.  We  first  quickly  show  well-known  results.  To  leading  order, 
y  =  yQ( v) ,  which  is  even  in  where 

[e* (T)]2y0  +vy(y0,T)=o  .  (1.2) 

♦'V 

Thus,  the  energy  E(T)  is  slowly  varying 


?[9*(T)]2(y0.  )2  +  V (y 0 , T )  =  E(T)  . 


(1.3) 


Here  v  is  the  fast  variable,  where  "j  We  will  eventually  determine 
the  phase  shift  o(T).  I.uke  [5]  showed  the  period  in  i \>  must  be  constant  and  hence 
the  frequency  ^(E,T)  satisfies 

.-1 


.(E,T)  =  9  1  (T  )  = 


,yo_. ..  % 


2  ]  'max 
%  .  S2  [E(T)_  V(yQ,T)] 


IT  i 


mi  n 


(1.4) 


where  yQ  (E,T)  and  y„  (E,T)  are,  respectively, the  minimum  and  maximum  of  the 
m  in  um  ax 

fast  nonlinear  oscillator.  After  performing  a  perturbation  expansion. 


y  =  y0  +  ey1  +  ....  we  obtain 


L(yi)  =  Ri  ,  i  =  1,2,3 


9  •  •  •  t 


(1.5) 


where  L  is  the  linearized  operator  resulting  from  (1.2)  and  where  is  the 


corresponding  right-hand  side  described  in  more  detail  in  §2.  By  the  Fredholm 


alternative  (for  (1.5)  with  i  =  l),  equivalent  to  eliminating  secular  terms,  we 
obtain  Kuzmak's  result  [4]  that  the  action  I  is  dissipated: 


dl 

df 


+ 


D  =  0 


(1.6) 


where  the  action  is  twice  the  average  value  of  the  kinetic  energy  divided  by 
the  frequency : 


I(E.T) 


CJ 


(E.T) 


r  1 
-0 


(y0  ) 

v 


(1.7) 


and  where  D  is  energy  dissipated  over  one  fast  oscillation: 


D  = 


fl 

•0 


(1.8) 


Equation  (1.4)  and  (1.6)  determine  E(T)  and  0'(T).  If  there  is  no  dissipation, 
then  the  action  I  is  conserved  (to  this  order)  and  is  the  well-known  adiabatic 
invariant. 

Next  in  §2,  we  determine  the  second-order  differential  equation  for  the 
phase  shift  c(T)  by  examining  the  0(e)  terms  of  the  exact  version,  developed 
by  Whitham  ([f]  and  [10]),  of  the  action  enuation.  Our  new  result  is 


d  I  do/dT  i 

di  l^rrrry  , 


=  0  . 


(1.9) 


Our  method  only  requires  a  calculation  of  y-|,  and  hence  is  much  simpler 

than  the  usual  elaborate  procedure  of  applying  the  secularity  condition  at 

0(c2).  We  relegate  this  to  Anoendix  A,  correcting  Dobrokhotov  and  Maslov  [1]. 

T n  (1.9),  u;E  is  the  rate  of  change  of  frequency  with  respect  to  energy  keeping 

T  fixed  via  (1.4).  Thus,  the  nerturbation  of  the  modulated  frequency  d$/dT 

divided  by  is  a  new  adiabatic  invariant.  Equation  (1.9)  easily  determines 

c (T)  from  initial  conditions  as  is  discussed  in  §3. 


Enuation  (1.9)  is  unexpected.  In  general,  perturbations  effect  the 
functional  form  of  the  dependence  of  the  frequency  on  the  amplitude.  However, 
in  §4  we  show  that  for  dissipative  perturbations  (limited  to  the  eh  term  in 
(1.1)  which  is  odd  in  dy/dt)  the  variation  of  both  frequency,  =  dq/dT, 
and  action  are  only  due  to  changes  in  the  averaged  perturbed  energy  E^ .  For 
example,  it  is  shown  that  io-|  =  E-j  ,  which  uses  the  leading-order  amplitude 

dependence  of  the  frequency.  Thus,  (1.9)  states  that  the  averaged  perturbed 
energy  does  not  vary.  Furthermore,  in  §4  it  is  shown  that  (1.9)  results  from 
considerations  of  the  linearization  of  the  equations  for  action  and  frequency 
with  resDect  to  changes  of  the  averaged  energy.  We  believe  the  simplicity 
of  (1.9)  is  due  to  damping  primarily  influencing  the  amplitude,  not  the  frequency 
However,  the  result  (1.9)  fundamentally  depends  on  the  damping  through  the 
dissipation  (1.6).  We  reiterate  that  action  is  not  an  adiabatic  invariant 
but  $'(T)/tj£  is.  Thus,  (1.9)  completes  the  description  of  the  leading  order 
problem  for  long  time;  it  shows  that  *(T)  is  not  necessarily  constant.  For 
other  tyoes  of  perturbations,  the  phase  shift  will  satisfy  nonhomogeneous 
versions  of  (1.9). 

2.  Stronqlv  Nonlinear  Oscillators. 

We  consider  the  second-order  nonlinear  differential  equation  which  describes 
strongly  nonlinear,  slowly  varying,  and  weakly  damped  oscillators: 


7  +  e  h  (y,  dt-*  T)  *  Vv(y,T)  •  0  , 
dr  01  y 


(2.1) 


where  V(y,T)  is  a  nonlinear  potential  depending  on  the  slow-time  scale  T  =et. 
We  assume  the  potential  admits  periodic  solutions  to  (2.1)  (for  all  fixed  T)  on 
the  fast  time  scale  t  without  dissipation  (i.e.,  when  e  =  0).  Also,  in  order  for 
the  small  perturbation  c  h  to  be  dissipative,  we  assume  h(y  ,  dy/dt,  T)  is  odd  in 
dy/dt.  This  allows  negative  dissipation  as  well  as  situations  such  as  the 
van  der  Pol  oscillator  in  which  both  positive  and  negative  dissipation  may  be 
present  simultaneously. 


We  use  the  method  of  multiple  scales  with  the  fast  time  scale  ip  and 
slow  time  scale  T: 


=  +4>(T) 


7  ■*  Et  , 


so  that,  for  example. 


{ft  -  [8'(T)*«*'<T)]$  ♦<£ 


In  the  analysis  to  follow,  we  will  rederive  0(T)  in  the  standard  way  and  correctly 
derive  the  phase  shift  <}>(T)  for  the  first  time.  In  this  manner,  (2.1)  becomes 
exactly 


(e*  +  ed>' )  +  e[(0"  +  e<t")y^  +  2(0'  +  e4>’)y^T]  +  eh 

+  c2yTT  +  Vy^y,T^  =  0  * 


(2.2) 


We  now  introduce  the  perturbation  expansion 


y(<KT)  =  y0(^,T)  +  ey-|  (>j>*T)  +  cZy2(^tT)  + 


It  is  convenient  to  note  that 


Vyty.T)  -  Vy(y0,T)  *  (cy,  ♦  t2y2)Vyy(y0,T)  +  e2  V(y0’T)  *  0(c3) 
h(y,v,T)  =  h(y,6'y^  ♦  ci'y^  +  eyT,T) 

=  h(yo,0'yo  .T)  +  e[h  v,  +  hv(<t>’y0  +  y«  +  0'y]  )]  +  0(e2)  , 

y  Hi  l  W 


where  v  e  dy/dt,  hy  =  hy(y0,6'y0  ,T)  and  hy  e  hv(yo,0'yo  ,T> 


The  result  of  substituting  the  perturbation  expansion  into  (2.2)  is  the 


nonlinear  oscillator  equation 


[e'(T)fy0  +  V  (y0,T)  =  0 

vv  y 


(2.3) 


and  its  higher  order  perturbations 


L(yk)  =  Rk  ,  k  =  1,2,3 . 


(2.4) 


where  L  =  (6 1 )  — *•  +  V  (yn,T)  is  the  linearization  of  the  nonlinear  operator  in 

dip  yy  u 

(2.3)  and 


-20Vy  -  0"y  -  20-yQ  -  h(yQ,e'yQ  ,T) 

%  %  UvT  0  % 


2e'4>y-|  -  d"y1  -  2e,y1  -  U>Ty0 

w  +>  Vr  w 


"  y0TT  ’  '  1  Vyyy^0,T^yl 

“  Vl  "  hvU'y0  +y0T  +  0,yl  }  * 

3  *  T  * 


We  now  examine  the  0(1)  equation,  (2.3).  Multiplying  by  yn  and  integrating 


once  yields  the  conservation  of  energy  equation 


\  (e')2(y0  )2  +  V(y0,T)  =  E(T)  , 


(2.5) 


•r 


uwimw  muwwtwm  wronw  m  m  vwnmm  tsxvwi 

K' 
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where  E(T)  is  slowly  varying  in  a  well-known  way  that  we  will  rederive.  Equation 

ay0 


(2.5)  defines  curves  in  the  (yg,  )  plane  for  fixed  T.  Because  of  our 


assumptions  concerning  the  potential,  these  curves  will  be  closed  and  yg 

will  be  periodic  in  ip,  oscillating  between  yQ  (E,T)  and  yQ  (E,T)  with 

umi n  umax 

V(yQ  ,T)  =  E  and  V(yg  ,T)  =  E.  Since  only  even  derivatives  with  respect  to 


min  max 

ip  occur  in  (2.3),  we  can  define  ip,  without  loss  of  generality,  such  that  yg(^,T) 


is  an  even  function  of  ip,  and  so  that  ip  =  0  corresponds  to  yg(0,T)  =  yQ 


and 


3yr 


min 


dip 


(0,T)  =  0  (see  Figure  1 ) . 


(E,T) 


Figure  1:  Phase  plane  solution  of  yn  for 


fixed  E,T. 


Using  quadrature  on  equation  (2.5),  we  obtain 


ip  =  0'  (T) 


ry0(^,T) 


dy  r 


Vn  [E-  V(y0,T2' 

mi  n  w 


The  period  of  the  oscillation  in  the  fast  variable  is  equal  to  the  loop  integral 
around  the  entire  curve,  or  equivalently: 


P(T)  =  29 ' (T) 


fyo_ 


dyf 


max 
■\ni  n 


yo„,_  v7  [E(T)  ■  y(y0’T)32 


(2.6) 


.Wii 


r-.i. 


i-M1 
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If  P  does  depend  on  T,  as  is  written,  then  as  Luke  [5]  showed 

•^■(Yq  U  +  nP(T),T))  =  n  Yq  (^>t)  +  (’KT),  which  is  unbounded  for  large  n. 

dP  ^  1" 

Thus,  jy  =  0  and  the  period  is  constant.  The  period  may  be  normalized  to  any 
convenient  constant,  and  we  will  choose  1  (although  2tt  is  often  used): 


j(E,T)  =  0  1  (T)  = 


cYq  (E,t) 


dyr 


,-l 


yo. 


max 

min 


(e.o  /2  [E(T>  -  v<Vr>]s 


(2.7) 


The  actual  period  of  the  oscillation  (to  first  order)  in  the  real  time  t  is 
1/e *(T)  SO  we  are  justified  in  referring  to  6 1 ( T )  as  the  frequency.  We  have 
introduced  the  notation  uj(E,T)  because  partial  derivatives  with  respect  to  E 
(keeping  T  fixed)  will  be  important.  (Note  the  difference  between  3w/3T 


and  dui/dT  =  3uj/3T  +  3^/3E  dE/dT.) 

In  order  to  examine  higher  order  equations,  we  first  need 
to  obtain  some  results  on  the  linearized  operator.  By  taking  the  derivative  of 
equation  (2.3)  with  respect  to  t j,  we  find  that  y^  is  a  homogeneous  solution  of 
L(u)  =  0.  As  was  stated  before,  y0(^,T)  is  an  even  function  of  ip,  which  means 

that  Yq  is  an  odd,  periodic  homogeneous  solution.  Since  Yg  is  periodic, 

o 

L(u)  =  0  is  equivalent  to  a  Hill's  equation.  From  Floquet  theory  the  form  of  the 
second  homogeneous  solution  is  -  A(T +  K(^,T)  where  K  is  even  and  periodic 

’P  gp  I 

in  47  with  period  one.  Luke  [5]  showed  that  A(T)  f  0  if  -g-p  f  0  in  equation  (2.6). 

Thus,  in  the  usual  case,  yQ  is  the  only  periodic,  homogeneous  solution.  In 

V 

order  for  solutions  to  (2.4)  to  be  periodic,  the  Fredholm  alternative  states  that 
the  right-hand  side  of  (2.4)  must  be  orthogonal  to  all  periodic  homogeneous 
solutions  of  (2.4).  This  gives  the  condition 

^Note  that  in  the  linear  case,  .-P/3E  =  0  for  all  E  and  there  are  two  periodic 
solutions:  cos  Zr<p  and  sin  Z~ .. . 
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y0  Rk  d*  =o  . 

Jo  \  K 


This  is  equivalent  to  eliminating  secular  terms  in  the  usual  way. 
For  k  =  1 ,  (2.8)  becomes 


[  2oi<t'y0  yQ  dip  +  2u  j  yQ  yQ  d^+-^  f  yQ  2  dip 
'O  V  Y  J0  V  <P  J0  Uip 


(2.8) 


+  h(y0»  w  yQ  ,T)y0  dip  =  0  . 

•'0  ip  <p 


The  first  integral  is  zero  since  it  is  the  integral  of  an  odd  periodic  function 
over  one  period.  The  next  two  integrals  may  be  combined  to  give: 


,y  [o)(E,T)  yn  dtj;  ]  +  h(yn,wyn  T)y~  d  ip  -  0 

dT  h  %  o  0  %  °<p 


(2.9a) 


This  is  a  well-known  result  [6].  We  define  the  action,  I(E,T),  to  be  twice 
the  average  value  of  the  leading  order  kinetic  energy  divided  by  the  leading 
order  frequency: 

I (E,T)  =  w(E,T)  [1(yQ  )2d*  =  2  fmax  ft  [E  -  V(yQ,T)f*  dyQ  ,  (2.10) 

®  v  ^mi n 

using  (2.5)  and  a  change  of  variables  (i.e.,  integrating  in  the  phase  plane). 

We  also  define  the  dissipation  (the  energy  dissipated  over  one  fast  oscillation): 


(2.11) 


D(E,T)  -  j  h ( y o » u, y g  » T )  y  g 
0  r 


Equation  (2.9)  becomes 


i 


E 


-  -  m  «  —  ,  w  —  -  «  ■>_  *■_  %  “»  ■%,  %  *%  .  %  "N*  .V  wV  j,'* 


WWW 


If  there  is  no  dissipation,  the  action  is  conserved  (to  this  order)  and  is  the 
well-known  adiabatic  invariant.  If  h  =  h(T)  ^  and  thus  h(yQ,u)y0  ,T)  =  h(T)wy 
then  (2.9b)  simplifies  to 


+  h(T)I  =  0  , 


-/'h(s)ds 

which  ([2], [4])  has  the  simple  solution  I(T)  =  l(0)e  .  In  any  case, 

however,  the  only  unknowns  in  (2.7)  and  (2.9b)  are  w(E,T)  and  E(T),  so  that 

these  equations  provide  a  closed  system  for  co  and  E. 

In  order  to  find  the  equation  for  (T) ,  it  will  be  necessary  to  examine 

higher-order  terms  and  thus  to  know  something  about  the  form  of  y^(ip,T).  By 

taking  the  partial  derivative  of  equation  (2.5)  with  respect  to  E,  noting  that 

^j(E,T),  we  obtain  L(yQ  )  =  -2uoj^y0  .  By  direct  substitution,  we  also  find 

E  W 

p 

that  L(’^yg  )  =  2  cj  yQ  .  Therefore,  yH  =  uyQ  +  is  a  homogeneous 

y  t p\l>  UE  y 

solution  to  (2.4),  which  is  even  in  Notice  that  yn  is  even  and  periodic, 

UE 

and  thus,  yH  is  in  the  form  predicted  by  Floquet  theory.  (Furthermore,  we 

note  that  there  are  two  periodic  solutions  if  =  0.  For  the  remainder  of  this 

paper  we  assume  (^(E.T)  f  0  since  we  believe  that  important  cases  correspond 

to  the  period  in  real  time  t  being  a  monotonic  function  of  the  energy.)  Using 

the  facts  that  yg  is  even  and  h(y,v,T)  is  odd  in  v  (which  means  h  is  odd  in 

allows  us  to  separate  R-,  =-2  :'yn  +  R.  where  R,  consists  only  of  terms 

1  'odd  'odd 


liA 


•TJV  V  V.^.VJV-.  «-.  « 
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=  .  1L 


-  "y°,  +  yin 

*  podd 


(2.12) 


where  a  specific  odd  part  of  the  particular  solution  may  be  obtained  by 
reduction  of  order 


yl  =  2 

podd  u 


'■P' 


Ri  .  d  ^ 


1^  Jo  ]odd  °5 


0  V,* 


^ —  d  p' 


We  will  not  need  this  expression  for  y^ 


Podd 


To  obtain  the  general  solution,  we  now  include  the  two  homogeneous  solutions, 
yQ  and  yH  given  above: 


y1  =  C1  (T)yQ  +  B-|  (T)  (uyQ  +  u ^yQ  )  -  ^  pyQ  +  y 


P  P 


odd 


In  order  for  y-j  to  be  bounded,  the  coefficient  of  ipy^  must  be  zero;  therefore, 
B^(T)  ■  <t'/uw£»  in  which  case 


yi  =  ci(T)  yo  *  *7%  *  yi„ 

P  t  p 


(2.13) 


odd 


1 


In  Appendix  A,  we  show  how  to  determine  an  equation  for  <t(T)  by 

p 

eliminating  the  secular  terms  in  the  '“'(e.  )  equation,  (2.8)  with  k=2.  Here, 
we  obtain  the  desired  result  in  a  simpler  way;  we  use  the  oeriodicity 
(secular)  condition  to  derive  an  action  eouation  valid  to  all  orders. 

Whitham  [10]  showed  that  an  exact  equation  for  action  may  be  obtained  by 

multiplying  (2.2)  by  y  and  integrating  over  using  the  assumption  that 

•¥ 

y  is  periodic  in  (with  neriod  1): 


4r  !  (•?'  +  e*')  '  y2d  .  +  f  h(y,v,T)y  d  y  =  0  . 

a]  -  Jo  -  J  JO  v 


(2.14) 


This  is  equivalent  tothe  usual  secularity  conditions  at  all  orders.  In 
particular,  to  leading  order  this  yields  (2.9).  More  importantly,  the 
phase  shift  will  be  determined  by  the  0(e)  terms  of  (2.14): 


[2.  ^ 

L  jo 


yo.  vi,  dy  +  *'  J0%  d 


fi  r  n  f1 

*  Jo  !V>  *  hv  *  y°T  +  9'V  %d(l  *  J  hy!4,  dt  '  0 


(2.15) 


Since  hy  is  even  in  c  and  hy  is  odd  in  \p,  only  the  even  part  of  y-j  is  needed 
in  (2.15).  From  (2.13) 


=  ill] 

'even  ~E 


(2.16) 


and  thus  (since  h^  y  T  yfl  d.  =  0  from  parity  considerations)  equation 


(2.15)  becomes 


+  I  y0  +  ^yn  )yn  h 


E 


E'O,  ,J0  "v 

V-  ^  V  V 


(2.17a) 


*  Vo  hy  +  yo  F  h  (J'o- "^0,  *T) 

L  U  WC.  U' 


d  -I  =0 


j 


d  r4> '  3  /  f1  2  .  »  -i  <J) 1  5 

ordT  ^  aE  UJ0y0/-)]  t;SE, 


h(y0,wy0  ,T)y0  d* 

0  ty  tp 


(2.17b) 


Equation  (2.17a)  is  equivalent  to  (A. 2),  derived  from  the  secular ity  condition 
2\ 


at  0(z  ).  Thus,  the  elaborate  calculations  in  Appendix  A  are  not  necessary 
to  obtain  the  phase  shift.  In  fact,  we  do  not  even  need  an  expression  for  R^. 
Using  the  definition  of  action  (2.10)  and  dissipation  (2.11),  we  obtain 


?r  'e]  *  dSI  de  ■ 0  • 


(2.17c) 


However,  action  is  dissipated  via  (2.9b),  D  =  -dl/dT,  so  that  (2.17c)  becomes 


d  r  $'(T)  t  "I  C '  (T)  5  /  d  I  \  _  « 

bt  arrfy  !e]  ■  ^TErfy  sr  W  - 0  • 


which  is  equivalent  to 
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r "ld*  dyn  , 

Since  from  (2.10)  IF  =  2  [E  -  V(yn,T)]  2  f  0  , 

L  J y  .  0  v7 


it  follows  that 


A  r  £0) 

dT  Lioe('E,1 


]  =  0  . 


(2.18) 


Thus,  (2.18)  defines  a  new  adiabatic  invariant.  Note  that  this  is  valid  for 


nonlinear  oscillators  with  rather  arbitrary  (smal 1)  dissipation.  An  alternative 

~^EE 

expression  occurs  if  we  observe  that,  from  above,  uv  =  .  Here  cuF  is 


determined  from  (2.7)  (do  not  differentiate  under  the  integral  sign),  where  we 


note  E(T)  has  been  determined  earlier  from  the  coupling  of  (2.7)  and  (2.9b). 


We  refer  to  <$(T)  as  the  phase  shift  and  $'(T)  as  the  perturbation  of  the 


modulated  frequency.  In  the  introduction  to  this  paper  we  have  described  the 


meaning  of  (2.18).  In  particular,  $(T)  and  <t>'(T)  are  now  easily  determined 


from  their  initial  conditions. 


Notice  that  the  only  unknown  in  (2.17)  or  (2.18)  is  <t>(T).  The  arbitrary 


coefficient  of  the  homogeneous  solution  in  (2.13),  C-j(T),  has  dropped  out  of  the 
equation  as  Luke  showed  [5].  We  note  that  it  is  possible  for  <J>'(T)  =  0,  as 
Kuzmak  [4]  assumed,  but  in  general  (2.18)  is  a  second-order  equation.  It  is 


important  to  recognize  this  for  the  arbitrary  initial  value  problem  (see  §3). 


We  can  determine  by  looking  at  the  0(e)  terms  of  the  action  equation  (2.14). 
Proceeding  in  this  manner,  we  can  find  Ck  by  looking  at  the  0(ck+l)  action 


equation.  It  is  very  imnortant  that  this  orthogonal i tv  condition  not  involve 


Cf.+^,  since  does,  in  general,  involve  yk+^.  Luke  [5]  showed  this  to 


true  for  the  non-dissipatlve  case  using  secularity  conditions,  and  it  is  easy  to 


show  for  the  general  case.  As  is  the  case  for  <+>(T),  we  would  find  that  Ck(T) 


WCWJ 


satisfies  a  second-order  equation. 


3.  Initial  Value  Problem 

The  importance  of  the  new  adiabatic  invariant  lies  in  the  relevance  of 
' (T )  (in  particular  <j>'(0))to  the  initial  value  problem  for  (2.1): 


2 

y ( 0 )  =  oiq  +  ca-j  +  e  o.£  +  ... 

^  (0)  =  60  +  e81  +  e262  +  ... 


(3.1) 


with  a.  and  g.  given.  By  using  the  multiple  scale  assumption  and  substituting 
the  asymptotic  expansion,  we  obtain 


since  ^(0)  =  $(0)  because  6(0)  =  0.  To  leading  order,  (3.2a)  for  i  =  0  and  (3.2b) 
form  a  closed  system  for  <b ( 0 )  and  E(0),  since  from  (2.7)  © '  ( T )  can  be  written 
as  a  function  of  E(T)  and  a  known  function  of  T  ( 6 ' ( 0 )  =  to(E(0) ,0) ).  E(T)  solves 
a  first  order  equation  (2.9a)  or  (2.9b),  and  therefore  (3.2a)  and  (3.2b)  determine 
E(T)  uniquely  (and  0'(T)  =  ..  (E,7)).  In  addition,  6 ( T )  is  uniquely  determined 
by  integration  since  9(0)  =  0. 


;'(T)  corresponding  to  the  leading  order  long  time  solution.  This  is  not  surprising 
as  oerturbed  initial  conditions  change  the  energy  and  hence  the  frequency. 


Using  the  expression  for  y^  ,  we  obtain  from  (3.3)  and  (2.18) 


odd 


'(T)  _ 
TT7TT  ' 


(o) 


u»evE(0),0)  =  1 


D 


yo  d 


-  “yo  yoT  + 

l  T 


'V 

hy0  d  * 
0  v 


+  U  I  e^Q  -  a1yQ  j 


(3.4) 


J 


T=0 
=  <t>  ( 0 ) 


To  simplify  this,  we  differentiate  the  energy  equation  (2.5)  with  respect  to  T 
and  inteqrate  it  from  v=0  to  j: 


i  dE  .  i  rv 

uj  CJ  J  Q 


VT  d  *  . 


(3.5) 


Evaluating  (3.5)  at  ,  =1,  yields  the  following  expression  for  dE/dT,  where 
the  action  equation  (2.9)  has  been  used: 


dE 

dT 


frhy0  d 


(3.6) 


By  substituting  (3.6)  in  (3.5)  and  by  noting  that  the  1 hs  of  (3.5)  appears  in 
(3.4),  we  obtain 


T'(T)  _  q'(0) 

■Ie(e,t)  ’  ~JT[oJ7o] 


*10)  f  f 

J0  1 


1  f  v. 


--  hyo 


+  UB^q 


Vo,  * 


<t> ( 0)  c  V 

—  -  hy 

0  *•  w 


0,1 


d  4,' 


T  =  0 

4>=4>(0)  . 


(3.7) 


Equation  (3.7)  is  significant.  It  determines  the  perturbed  frequency  and  shows  how 
it  depends  on  the  dissipation,  Vvq  d  C’,  the  slow  time-dependence  ofthe  potential,  Vy, 
and  the  perturbed  initial  conditions,  and  6-| .  For  example,  a  simpler  expression 
results  if  =  ?^  =  0  and  there  is  no  dissipation  (h=0);  the  frequency  changes  are 
only  due  to  a  time-deoendent  potential.  If  in  addition,  Vy  =  0,  then  q(T)  is  constant, 
as  is  clear  since  in  this  snecial  case  there  is  no  slow  variation  of  the  solution. 
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4.  Discussion:  Higher-order  Frequency  and  Averaged  Perturbed  Energy. 

The  leading-order  equations  for  the  slow  variation  of  a  nonlinear  oscillator 
are  the  amplitude-dependent  frequency  (2.7)  (based  on  the  leading-order  energy) 
and  the  action  equation  (2.9).  Our  result  concerning  the  phase  shift  is  best 
understood  by  considering  higher-order  expressions  for  both.  In  this  section, 
we  note  that  changes  in  the  frequency  and  action  equations  only  result  from 
the  linearization  of  them  due  to  perturbing  the  averaged  energy,  if  the 
perturbation  is  dissipative,  as  we  have  assumed. 

In  Appendix  B,  it  is  shown  that  the  frequency  (2.7)  is  valid  to  0(e) 
if  the  average  perturbed  energy  E  is  utilized  in  (2.7).  Thus,  the  Taylor 
expansion  of  the  frequency  u  =  ^(E)  yields 

t-’l  ~  E  i  »  ( 4, 1 ) 

where  E-j  is  the  0(e)  perturbation  of  the  averaged  energy  and  to,  =  d^/dT  is 
the  0(e)  perturbation  of  the  frequency  (the  frequency  shift  or  derivative  of 
the  phase  shift).  Here,  frequency  changes  are  only  due  to  the  averaged 
energy  changes,  and  thus,  we  might  reason  intuitively  that  (4.1)  should  be 
valid. 

Equation  (4.1)  is  quite  significant.  For  example,  the  leading-order 
nonlinear  oscillator  satisfies 

yQ  =  yQ(-:;  E(T),  T)  . 

We  expect  that  some  of  the  perturbed  solution  y-|  results  from  changes  in  E. 

In  fact. 


& 
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*1 

even 


Ei  yoc 


(4.2) 


since  from  (4.1)  E-j  =  d0/dI/cuE ,  showing  that  (4.2)  is  equivalent  to  (2.16). 
Thus,  (2.16)  or  (4.2)  corresponds  to  the  response  due  to  perturbed  average 
energy  levels. 

Furthermore,  (4.1)  leads  to  an  understanding  of  the  phase  shift.  Since 
action  is  valid  to  higher-order  from  (2.14)  and  both  action  and  the  dissipation 
depend  on  the  averaged  energy  E  [7],  it  follows  that 


^  I (E,T)  +  D(E,T)  =  0 


(4.3) 


To  leading  order,  this  yields  (2.9).  Higher-order  terms  may  be  simply 
obtained  by  perturbing  the  averaged  energy: 

^  (E1  IE)  +  E1  De  =  0  •  (4.4) 

However,  the  partial  derivative  of  (4.3)  with  respect  to  E  is  also  valid. 
Thus,  (4.4)  becomes 


I 


_d_ 
E  dT 


0  . 


This  shows  that  the  averaged  perturbed  energy  E-|  does  not  vary.  Using  (4.1), 

we  obtain  the  adiabatic  invariant,  —  d^/dT: 

WE 


d  (d<£/dT j 


0  . 


(4.5) 
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Our  results  are  only  valid  for  the  dissipative  perturbation  ch  .  For  non- 

p 

dissipative  perturbations  or  0(e  )  dissipative  perturbations,  (4.5)  becomes 
a  non-homogeneous  equation  for  d$/dT. 

Recently,  we  have  extended  the  ideas  of  this  paper  to  obtain  the  modulated 
phase  shift  for  oscillatory  solutions  of  nonlinear  dispersive  waves 
([7], [8],  and  [9]). 
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In  this  appendix,  we  show  that  the  higher-order  action  equation  (2.15) 
is  eouivalent  to  the  result  of  eliminating  the  secular  terms  in  the  0{c2)  equations 


R2y  d^  =  0 

0  ip 


(A.l) 


fl  2 

The  term  I  vvvv(yn»T)yi  y0  in  (A.l)  appears  to  be  difficult  to  calculate 

jo  yyy  u 

since  y-|  is  given  by  (2.13).  However,  if  we  integrate  by  parts,  recall  the 
definition  of  L,  and  then  use  L(y-j)  =  R-j,  we  obtain 


r 1  2 

j  „V>(,0JI,1  yo.d4, 

u  v 


=  ^LVV^i  ^ 


*  -  2  CKy^  -  to  y1  ]y,  d  ip 

■>0  yip  \p 

f1 

=  +  2  C2  to  <t  +gfy0  +2wy0  +h(y0»^0  »T)  lv1  d, 

J0  VV  'P  ^T  ij; 


where  the  periodicity  of  many  expressions  has  also  been  used.  In  this  manner, 
from  (A.l),  we  obtai n 


20 L  ■^(yOyl,)d^+  2dTtu  LyO,yl,d*] 


0  "p  'ip 


fl  fl  2  fl 

+  h(y0.<-y0  ’T^yl  d'^  +  ^0  '  hvd^  +  w  yl  y0  hvd^ 

J0  v\p  \p  •’0  ip  -*0  ^  i> 

+  f  yn  yi  hu  do  +  0’  f  (yn  )2^3 
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This  corrects  an  error  in  Dobrokhotov  and  Maslov  [1]  due  to  an  omitted  term. 
The  first  integral  is  exact  and  vanishes  due  to  periodicity  considerations. 
The  underlined  integrals  vanish  from  parity  considerations.  Thus,  (A. 2)  is 
equivalent  to  (2.15)  derived  from  the  higher-order  action  equation. 


Appendix  B:  Averaged  Perturbed  Energy. 

It  is  useful  in  interpreting  the  phase  shift  to  show  that  (4.1)  is 
valid,  that  is,  to  evaluate  the  averaged  perturbed  energy.  The  averaged 
energy  is  defined  to  be 


E 


cH 


(B.l) 


When  y  is  expanded  and  the  multiple  scale  assumption  is  used,  the  leading- 
order  satisfies  (2.5)  with  E  replaced  by  Eg.  If  is  defined  as  the  0(c) 
perturbation  of  the  averaged  energy,  then  from  (B.l)  we  obtain 


f1  -  2 

e1  =  !  !e,%:(?'yi;  +  c ,yn.  +  ynr  -  (?')  y,y 


% 


'OT 


r  OiK* 


] 


(B.2) 


using  (2.3).  We  note  that  only  y,  is  needed  in  (B.2),  and  thus,  using 

even 

parity  arguments  and  (2.16),  we  obtain 


E1  =  ds/dT/.E  , 

proving  that  (4.1)  is  valid.  In  a  similar  way  it  can  be  shown  that  (4.3)  is 
valid  to  one  higher  order. 


References 


[1]  S.  Yu.  Dobrokhotov  and  V.  P.  Maslov,  Finite-zone,  almost-periodic  solutions 
in  WKB  approximations,  I  tog i  Nauki  i  Tekhniki,  Seriya  Sovremennye  Problemy 
Matematiki,  15  (1980),  pp.  3-94.  J.  Sov.  Math.  15  (1980),  1433-1487. 

[2]  J.  Kevorkian  and  J.  D.  Cole,  Perturbation  Methods  in  Applied  Mathematics, 
Springer-Verlag,  New  York,  1981. 

[3]  S.  Kogelman  and  J.  B.  Keller,  Asymptotic  theory  of  nonlinear  wave  propagation 
SIAM  J.  Appl .  Math.  24  (1973),  pp.  352-361. 

[4]  G.  E.  Kuzmak,  Asymptotic  solutions  of  nonlinear  second  order  differential 
equations  with  variable  coefficients,  Prikl.  Mat.  Meh.  23  (1959),  pp.  515-526 
J.  Appl.  Math.  Mech.  23  (195"),  pp.  730-744. 

[5]  J.  C.  Luke,  A  perturbation  method  for  nonlinear  dispersive  wave  problems, 
Proc.  Roy.  Soc.  London  Ser.  A  292  (1966),  pp.  403-412. 

[6]  G.  B.  Whitham,  Linear  and  Nonlinear  Waves,  John  Wiley  &  Sons,  New  York,  1974. 

[7]  R.  Haberman  and  F.  J.  Bourland,  Variation  of  wave  action:  Modulations  of 
the  phase  shift  for  strongly  nonlinear  dispersive  waves  with  weak  dissipation 
submitted. 

[8]  F.  J.  Bourland  and  R.  Haberman,  Dissipatively  modulated  nonlinear  dispersive 
waves  with  a  Lagrangian  formulation,  in  preparation. 

[9]  R.  Haberman,  The  modulated  phase  shift  for  weakly  dissipated  nonlinear 
oscillatory  waves  of  the  Korteweg-deVries  type,  submitted. 

[10]  G.  B.  Whitham,  Two-timing,  variational  principles  and  waves,  J.  Fluid  Mech. 
44(1970),  pp.  373-395. 


i* | 


UPOfiR.TR*  87"  156  9 


The  Modulated  Phase  Shift 


for  Weakly  Dissipated  Nonlinear  Oscillatory 
Waves  of  the  Korteweg-de  Vries  Type 


Richard  Haberman 


Department  of  Mathematics 
Southern  Methodist  University 
Dallas,  Texas  75275 


Revised  September  1987 


Abstract 


Nonlinear  dispersive  oscillatory  waves  are  analyzed  for  Korteweg-deVries 
type  partial  differential  equations  with  slowly  varying  coefficients  and  arbitrary 
small  perturbations.  Spatial  and  temporal  evolution  of  the  amplitude  parameters 
are  determined  in  the  usual  way  by  the  possible  dissipation  of  the  wave  actions 
for  both  momentum  and  energy.  For  dissipative  perturbations,  both  wave  actions 
are  shown  to  be  valid  to  a  higher  order.  Thus,  the  first  variation  of  the  wave 
action  equations  is  used  to  derive  equations  for  the  slow  modulations  of  the 
phase  shift.  It  is  shown  that  the  phase  shift  satisfies  a  universal  set  of  two 
coupled  equations,  each  independent  of  the  small  dissipative  perturbation. 


Introduction 

Kuzmak  [4]  made  fundamental  contributions  to  the  method  of  multiple  scales 
in  his  analysis  of  the  effects  of  weak  dissipation  and  variable  media  on  strongly 
nonlinear  oscillators  satisfying  ordinary  differential  equations.  The  slowly 
varying  frequency  was  shown  to  equal  the  local  frequency,  depending  on  the 
amplitude.  Kuzmak  [4]  also  obtained  the  differential  equation  for  the  slowly 
varying  amplitude  (a  nonlinear  generalization  of  the  dissipation  of  action).  For 
strongly  nonlinear  dispersive  waves  (but  without  dissipation),  satisfying  the 
Klein-Gordon  equation  (a  partial  differential  equation)  with  variable  coefficients, 
Luke  [5]  obtained  the  dispersion  relation  and  the  partial  differential  equation 
for  the  spatial  and  temporal  slow  evolution  of  the  amplitude.  Whitham  [6] 
generalized  this  by  intr^Jucing  the  concept  of  wave  action,  enabling  as  an  example 
the  calculation  of  the  effects  of  variable  media  and  weak  dissipation  on  the 
cnoidal  waves  for  the  Korteweg-deVries  equation. 

Kuzmak  [4]  assumed  the  phase  shift  was  constant,  while  for  a  dissipationless 
system,  Luke  [5]  claimed  that  the  phase  shift  could  be  constant,  an  important 
distinction  that  has  been  very  understandably  overlooked  by  many.  Unfortunately, 
when  Dobrokhotov  and  Maslov  [2]  reinvestigated  this  issue  they  made  a  critical 
error.  Although  they  obtained  a  second-order  equation, 

their  results  were  incorrect.  Furthermore,  they  only 
analyzed  unperturbed  dissipationless  systems. 

Quite  recently,  Bourland  and  Haberman  [l]  analyzed  nonlinear  oscillators 
with  dissipation.  Using  some  of  the  good  ideas  of  Dobrokhotov  and  Maslov  [2], 
a  simple  second-order  ordinary  differential  equation  for  the  phase  shift  was 
obtained.  Most  significantly,  it  was  of  universal  form,  independent  of  the  type 
of  dissipation.  Furthermore,  there  was  a  simple  adiabatic  invariant,  enabling 
the  phase  shift  to  be  determined  for  arbitrary  linear  or  nonlinear  dissipation. 


|Ly.  **.  •  y.  •*.  ..‘  V 
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The  simplicity  and  generality  of  the  phase  shift  for  nonlinear  oscillatory 
ordinary  differential  equations  strongly  suggested  that  similar  results  would 
be  valid  for  weakly  dissipated  nonlinear  dispersive  waves  propagating  in  a 
slowly  varying  media.  Efforts  were  simultaneously  undertaken  for  Klein-Gordon 
and  Korteweg-deVries  type  equations;  each  guided  the  development  of  the  other. 
For  the  Klein-Gordon  equation,  the  dissipation  of  wave  action  is  fundamental, 
in  one  spatial  dimension. 


1*  +  =  _r 

IT  3X  ’ 


(1.1) 


where  I  is  the  wave  action,  q  is  the  flux  of  wave  action,  and  R  is  the 
dissipation.  Haberman  and  Bourland  [3]  showed  that  the  phase  shift  satisfies 
the  same  partial  differential  equation,  independent  of  the  type  of  dissipation: 


D-j-I  +  Dxq  =  0  , 


(1.2) 


where  D  is  the  Taylor  series  operator 

resulting  from  changes  in  the  wave  number  and  averaged  perturbed  energy. 

A 

In  [3]  it  is  shown  that  D  depends  on  the  phase  shift  in  a  linear  way,  so  that 
(1.2)  is  a  linear  partial  differential  equation  for  the  phase  shift.  Once 

^  A. 

D  has  been  determined,  (1.2)  can  be  obtained  from  (1.1)  by  replacing  by  D-j- 

,  A 

and  ~  by  Dx  and  by  neglecting  the  dissipation. 

In  this  paper,  we  analyze  Korteweg-deVries  type  equations  for  a  variety  of 
reasons.  It  is  an  important  type  of  equation  because  of  its  frequent  occurrence 
in  physical  problems  and  it  can  have  great  mathematical  interest  in  the 
remarkable  cases  where  it  is  integrable  and  known  to  have  exact  multi-phase 
periodic  solutions.  However,  we  investigate  slow  variations  because  nonlinear 
periodic  waves  satisfy  a  third-order  equation  (not  self-adjoint)  which  is 


thus  more  difficult  than  the  corresponding  problem  for  the  Klein-Gordon  equation. 

For  Korteweg-deVries  type  equations,  there  are  two  independent  amplitude  parameters, 
while  for  the  Klein-Gordon  there  was  only  one. 


In  §2,  the  method  of  multiple  scales  is  used  to  analyze  Korteweg-deVries 
type  eauations  with  arbitrary  small  perturbations.  Two  equations  are  derived 
representing  the  well-know  dissipation  of  wave  action  for  momentum  and  energy, 
both  in  the  form  of  (1.1),  with  I,  q  replaced  by  I.,  q.  for  i  »  1,2.  In  §3, 
higher  order  terms  are  shown  to  be  obtained  from  the  first  variation  of  a  form  of 
the  wave  action  equations,  an  apparently  general  result  since  it  is  also  valid 

A 

for  the  Klein-Gordon  equation  [3].  In  §4,  one  operator  D,  depending  on  the 
phase  shift,  is  determined  such  that  (1.2)  is  valid  for  both  wave  actions.  In 
this  way,  the  phase  shift  for  Korteweg-deVries  type  equations  is  shown  to  satisfy 
a  coupled  system  of  partial  differential  equations  (1.2)  of  universal  form, 
independent  of  the  type  of  dissipation.  The  equations  (1.2)  are  easy  to  express 
once  the  operator  D  is  determined.  In  §5,  we  show  that  these  results  correspond 
to  the  linearization  of  the  dispersion  relation  and  the  wave  action  equations  with 
respect  to  perturbations  of  the  wave  number  and  the  averaged  values  of  the  amplitude 

/V 

parameters.  The  operator  D  is  now  the  Taylor  series  operator  in  these  parameters. 

2.  Korteweg-deVries  Type  Equations. 

In  this  paper,  we  analyze  nonlinear  partial  differential  equations  with 
slowly  varying  coefficients  of  the  Korteweg-deVries  (KdV)  type: 


u.  +  i  Q  (u,X,T)u  +  a^(X,T)u 
t  2  yuu'  ’  '  x  x  '  xxx 


=  - £  h(u,ux,u<A ,X,T)  -  e  hj (u,ux,uxx,X,T)  , 


(2.1a) 


where  ^  Guu(u»x>"0  is  the  convective  velocity.  To  correspond  to  nonlinear  wave 
phenomena  in  a  slowly  varying  media,  we  allow  quantities  to  depend  on  the  slow  time 
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and  spatial  scales,  T  =  ct  and  X  =  ex.  As  we  will  show  later  in  this  section, 

it  is  possible  to  assume  that  ^  Quu(u,X,T)  can  be  chosen  such  that  there  are  nonlinear 

oscillatory  wave  solutions  of  (2.1a)  (for  all  fixed  X  and  T)  on  the  fast  scale 

2 

t  when  the  perturbations  eh  and  e  h-j  are  absent.  Typical  examples  are 
1  3 

Q  =  -g-8  (X,T)u  ,  in  which  case  (2.1)  becomes  the  slowly  varying  and  perturbed 

Korteweg-deVries  equation,  and  Q  =  ^-8  (X.T)u^,  in  which  in  a  similar  way  (2.1a) 

relates  to  the  modified  Korteweg-deVries  equation. 

2  2 

The  term  e  h^  is  introduced  since  0(e  )  perturbations  can  influence  the 

leading  order  solution  for  long  time.  At  first,  we  calculate  the  effects  of 

2 

arbitrary  slowly  varying  perturbations  ch  and  e  h-j.  However,  our  primary  interest 

is  in  dissipative  perturbations.  Thus,  later  in  the  paper  we  assume  eh  is 

dissipative,  but  simplifications  will  only  occur  if  we  assume  in  addition  that 
2 

e  h|  is  non-dissipative. 

We  briefly  review  some  elementary  ideas  concerning  conservation  laws.  We 
interpret  u  as  momentum-density,  so  that  (2.1a) itself  expresses  the  possible 
dissipation  of  momentum  by  the  perturbations: 


It  *  &  <X  +  Ac*1  =  \  Qux +  2c“Vxx  - Eh  -  • 


(2.1b) 


1  2 

where  -jQ  +  a  u  is  the  momentum  -flux.  Here,  momentum  also  may  change  due  to 

c  U  A  A 

the  variable  media,  represented  by  the  terms  -iq  and  2aa  u  .  Similarly, 

Cm  HA  A  AA 

multiplying  (2.1a)  ty  u  shows  how  the  energy  may  be  dissipated: 


£<!u2>  *  + 


3x  3x 


3  2t 
7ux] 


-  chu  -  c  h^u 


luQux  ■  ?Qx 


(2.1c) 
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Thus,  iuQu  -  +  :s2[-^-  (uux)  -  |u2  J  is  the  energy  flux. 

Typical  linear  perturbations  are  h  =  hgU,  which  dissipates  both  momentum 

and  energy,  and  -h  =  vuxx,  a  diffusive  term  which  conserves  momentum,  but 

dissipates  energy.  In  this  latter  case,  (2.1a) is  of  the  KdV-Burgers  type. 

However,  for  the  moment  we  will  be  quite  general  and  make  no  assumptions 

2 

concerning  ch  or  e  h-j.  Later,  when  we  assume  that  eh  represents  dissipation, 
we  will  allow  any  kind  including  nonlinear  types. 

To  obtain  the  long  time  behavior  of  nonlinear  dispersive  waves  including 
the  ffects  of  a  slowly  varying  media  and  the  small  perturbations,  the  method 
of  multiple  scales  will  be  used.  We  follow  Luke's  [5]  procedure  for  the  Klein- 
Gordon  equation  as  re-examined  by  Haberman  and  Bourland  [3]: 


+  *(x,t; 


T  =  et 


X  =  cx 


This  method  is  particularly  suited  for  determining  the  modulations  of  the 
phase  shift  $(X,T).  For  example, 


h  =  (ex  +  e*xK  +  CU) 


In  this  way,  (2. la)  becomes  "exactly" 
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<eT  +  ‘V'V  +  gut  *  7<ex  +  £»x>  i  Du  *  qu  -  f  Q 


uX 


*  c  (X,T)[(6X  +  c»x)  Vk  +  3e(9X  +  E*x>  U»X 


*  5e(9xx  *  ^xx)(ex  +  E«x)  ♦  3e2  (9x  +  E*x)  uXXJj 


+  3e2(9xx  +  e*xx)ua  +  e2(Sxxx  +  E4.m)u^  *  e3uxxx] 


■ch[u,9  u  +  c$vu  +  eu¥ ,  (eY)u  ,+ 

A  y  A  ]{'  X  X  (j/lj/ 


+  2cVxu^  +  2c6x%x  +  £0xx%  +  0(e  X*T1 


- 2h -j [u ,  -xu  +  0(e),  (6x)2u  +  0(e), X.T]  . 


The  substitution  of  the  perturbation  expansion, 


u  =  uo  +  'ul  +  •  u2  +  •••  • 


will  be  facilitated  by  noting 


(2.2) 


Qu(u,X,T)  =  Qu(u0,X,T)  +  (£U]  +  c  u2)Qu(j  (u0,X,T) 
2 


2  U1 


♦  'VO  f  un,X,T)  ♦  0(  ‘  J) 
2  uuu  0 
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h[as  in  (2.2)]  =  h(Ug,9xUg  ,(9X)  uQ  ,  X,T) 


w 


+  e[huUl  +  hu  (exul  +  Vo  + 

X  r  A 


+  hu  ^Sxul  ,  +  2VxU0  +  29XU0  y  +  6XXU0  ^  +  °(£  )  » 

xx  ipi^/  i^X  ^ 


where  h  ,  ,  and  hu  are  evaluated  at  Ug,  6xUg  ,  Ug  ,  X,T. 

X  xx  \p  Vy 

We  obtain  in  this  way,  to  leadina  order,  the  nonlinear  ordinary  differential 

equation  which  represents  travelling  wave  solutions  of  the  KdV  type  equations: 


~TUG  +  2  “X  i,  VU0,X,T^  +  J  u0  =  0  ’ 


(2.3) 


2  2 

where  *  i  (X,T).  Equation  (2.3)  is  solved  in  the  usual  way.  Integration 
yields 


Tu0  ‘  KW’T>  *  'S>\  *  -B  '  (2'4) 

¥  i‘ 

where  B  B(X,T).  By  multiplying  by  2Ug  and  integrating,  we  obtain  an  equation 
which  can  be  interpreted  as  conservation  of  energy  for  the  system  (2.4): 

Vo  *  V(VX,T)  +  l2(“x)3(uo  )2  =  "  2Bu0  +  2A  ’  (2,5) 

V 

where  A  =  A(X,T).  The  parameters  A(X,T)  and  B(x,T)  are  similar  in  spirit  to 
those  introduced  by  Whitham  [6]  for  the  slowly  varying  KdV  equation.  Our 
analysis  is  valid  for  those  Q(un,<,T)  for  which  there  are  periodic  solutions,  see 
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2 

the  "potential"  5yUg  +  SxQ(ug,X,T)  of  (2.5).  From  (2.5),  we  can  define 
Uq(v.X.T)  to  be  an  even  function  of  We  introduce  the  wave  number  k  e  and 
frequency  *  =  -  ~ j  .  Luke  [5]  showed  that  the  period  in  the  fast  variable  y 
is  constant.  If  we  normalize  that  period  to  1 ,  we  obtain  the  amplitude 
dependent  dispersion  relation,  u  =  u(k(X,T),  A(X,T),  B(X,T),  X,T),  for  the 
nonlinear  dispersive  wave  in  a  nonhomogeneous  media: 


mi  n 


(2.6) 

where  the  wave  oscillates  from  Uq  to  Uq  ,  two  successive  zeroes  of  the 

min  max 

denominator  of  (2.6).  Here,  there  are  two  amplitude  parameters  A(X,T)  and 
B ( X , T )  unlike  the  case  of  the  Klein-Gordon  equation  [3]  where  there  is  only 
one  E(X,T) . 

The  higher  order  equations  in  the  perturbation  expansion  are 


L(ui)  =  Ri  i  =  1,2,... 


(2.7) 


where  L  is  the  linearization  of  the  nonlinear  operator  in  (2.3): 


1  O  0  ~  J 

L  -  ^t  +  rx-rQuu  +  ^x77 

ay 


(2.8) 


In  a  straight-forward  manner,  we  obtain  ,  the  term  we  will  need 

to  determine  the  amplitude  parameters  A(X,T)  and  B(X,T)  and  the  phase  shift 
:(x,T).  We  do  not  need  P-,  though  it  will  be  used  in  Appendix  A  to  show  that 
our  procedure  is  equivalent  to  t^at  jsed  by  Dobrokhotov  and  Maslov  [2]: 


9 


R1  ’  "Vo,  '  uOT  "  2  *x  dl  Qu  ‘  2  3X  Qu  +  2  QuX 


a  (3^x  <t)x  u0  +  35XU0.  +  3exX0XUO,|) 

W  \p  \pipX  y  \p 


“  ^(Uq*^^Uq  Uq  «  X,T) 


(2.9) 


,  =  -0T-U-.  -  u,  -  iev^-  (u?  Q  ) 

2  T  I  .  IT  2  X  3^  '  i  wuuir 
V  I 


‘  2  ^X  3y  ^lQuu^  ‘  2  3X  ^u1Quu^  +  2  U1  QuuX 


-  a2(36xixui  +39x<J)Xu0.i.  ,  +  39XU1  . 

V  y  y  V  y ’y  ^  y  A 


+  6Vxuo  .  +  3-xx“xui  +  3*XX6Xu0 


+  3“xx;,x“o  *  3exuoy  +  3exxuo.  *  6XXXu0  * 

f'  y  A  A  y  y  A  y 


hl  ^ u0’  jxu0  *GXU0  ,  ,X,T^  "  Ulhu  "  ^xul  +  *XU0.  +  UOv)hu 

y  ly  V  y  y  a  x 


2 

(ev  u,  +  2rvi>vun  +  2ev  un  +  ~  un  )  h 
XI  X  X  0  X  0  v  XX  0  '  u 

yy  y X  v  xx 

(2.10) 


Equation  (2.7)  is  a  third-order  equation.  In  order  for  u.  to  be  periodic 
(in  y),  the  Fredholm  alternative  implies  that  must  be  orthogonal  to  all 
periodic  solutions,  1  and  un,  of  the  homogeneous  adjoint  equation: 


Ri  d; =  0  i  =  1 ,2,  ... 


(2.11) 


Ri  uQ  &'j  -  0  i  =  1 ,2,  . . 


(2.12) 


At  0(e),  these  yield  the  well-known  equations  to  determine  A(.X,T)  and 
B( X ,T ) .  Using  (2.9)  and  (2.11),  di rectly  yields 


vr  [uod«  +  2  i*  {,V*  =  2  j0Quxd’> 


h(uo,0xuo  »X,T)di^ 

0  \p  \pip 

(2.13) 


since  is  periodic.  This  shows  that  wave-momentum  action  may  not  be 
conserved  due  to  the  rhs  of  (2.13).  Using  (2.9)  and  (2.12),  after  some 
elementary  algebra,  we  derive  a  similar  result  for  the  possible  variation  of 
wave-energy  action: 


•T  J0  2  u0  ij  *  2  ;x  j0  <u0°u  '  *  +  2  ■  2  Q  uQuX  d  ^ 


■2-t|^X>2  J0(uo/,J*] 


;  h(u0,  0  u0  ,  6X  u  ,X,T)u0d»  . 

-0  'i>  w 


(2.14) 


Here,  0^  is  the  partial  derivative  of  Q  with  respect  to  X  keeping  u  fixed. 

In  many  problems,  the  convective  velocity  Quu  only  depends  on  u,  and  thus  in 
this  case,  Qx  =  0,  preserving  the  conservation  form  of  (2,13)  and  (2.14)  if  h  =  0 
The  integrals  in  (2.13)  and  (2.14)  depend  on  A(X,T),  B(X,T),  k(X,T),  and 
- ;A,6,k,*,T) .  Thus  the  dispersion  relation  (2.6)  and  conservation  of  waves. 


-11- 

■^k  +  ~  w  =  0  ,  (2.15) 

combined  with  (2.13)  and  (2.14)  determine  the  amplitude  parameters  A(X,T) 
and  B(X,T)  and  the  dominant  phase  parameters  k(X,T)  and  w(A,B,k,X,T) ,  but  not 

the  phase  shift  <j>(X,T).  Equations  (2.13)  and  (2.14)  agree  with  Whitham's  [6] 

1  3 

result  for  the  KdV.  However,  for  the  KdV,  Q  =  ■jBu  and  thus  the  momentum 

2 

flux  Qu  is  proportional  to  the  energy  density  u  ,  reducing  by  one  the  number 
of  independent  integrals  in  (2.13)  and  (2.14).  These  results  also  show  that  even 
though  the  KDV  has  an  infinite  number  of  conservaiton  laws,  momentum  and  energy 
are  the  only  appropriate  ones  to  consider  when  perturbations  are  applied. 


3.  Higher-Order  Wave  Actions. 

In  Appendix  A,  we  determine  the  equations  for  the  modulations  of  the 

o 

phase  shift  c(X,T)  by  eliminating  secular  terms  at  0{t  ): 


fl 

j  R«  d  0  =  0  and 

J0 


r  1 

RpUndvp  =  0  . 

J0  2  0 


( 3 . 1  a ,  b ) 


This  approach  is  somewhat  involved  since  it  requires  R2  given  by  (2.10). 

Here,  we  use  a  simpler  method  which  yields  equivalent  results.  We  show  that 
the  phase  shift  is  determined  from  the  0(e)  wave  action  equations. 

As  Whitham  [6]  has  shown  for  Klein-Gordon  equations,  exact  equations  for 


the  wave  actions  may  be  derived.  For  KdV  type  equations,  an  exact  expression 
for  wave-momentum  action  is  obtained  by  integrating  (2.2)  over  \p  from  0  to  1 
and  applying  the  periodicity  of  u: 
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fl  i  3  (1  i  ft  9  ? 

3l  j„ud*  +  2  3X  [V*  -  2  ^uX  d  ^  +  E  “  ^3  Ju"* 


2  3X  JQ  ’u 


h  d  .  -  e 


J0 


h-,dr 


(3.2) 


By  first  multiplying  (2.2)  by  u  and  then  integrating,  we  obtain  (with  more 
effort)  an  exact  expression  for  wave-energy  action: 


TT  |Q  \  y2  d  *  +  7  Ik  [  (uQu  "  Q)dv  +  1  j0  Qxd'^  ’  2  |0  uQuXd^ 

-  I  a2  ^7  [(ev  +  cfv)2  f  u2d”]  -  3a2  [(ex  +  €<PX)  [  u  uxd^] 


9X  L v  X  +Xy 

fl 

Jo 


■'0 


2  2  f 

+  £  a  j  uuxxx  d  y  =  -  j  hudv  -  e  j  h^  u  d'^ 


(3.3) 


It  is  perhaps  even  simpler  to  obtain  (3.2)  and  (3.3)  by  directly  introducing 
the  multiple-scale  assumption  into  (2.1b)  and  (2.1c),  and  then  integrating. 
Wave-momentum  action  (3.2)  is  "identical"  to  the  leading-order  result  (2.13) 
if  the  entire  u  is  used  as  well  as  the  total  wave  number,  ex  +  e<}>x,  except 
for  the  effects  of  the  perturbations  and  0(e  )  dispersive  term.  However, 
wave-energy  action  (3.3)  is  alterred  slightly  by  the  additional  term 


r  1 


-  3a  C(ex  + ed>x)  u,  uxd  v]  .  This  is  an  0(e)  correction,  but  since 
Ug  is  even  in  it  is  in  fact  smaller  being  at  least  0(e  ).  Thus,  we  have 

shown  that  both  wave  actions  are  valid  to  at  least  one  higher  order  if  the 

2 

effects  of  the  perturbations,  eh  and  e  h^  are  properly  introduced. 


V ft.*' O’ 


iL~lO«  lib  ■>>  A*.  >««  a»|  |J  L*  L‘ 


These  exact  wave  action  equations,  (3.2)  and  (3.3),  are  equivalent  to 
eliminating  secular  terms  at  all  orders.  For  example,  (3.2)  is  the  same  to 
leading  order  as  (2.13),  and  (3.3)  relates  similarly  to  (2.14).  Moreover, 

evaluating  the  exact  action  equations  to  0(e)  is  equivalent  to  eliminating 

2 

secular  terms  at  0(e  )  (see  Appendix  A).  In  particular,  the  0(e)  wave- 
momentum  action  terms  from  (3.2)  yield 


3  f1  i  .  f1  ]  rl  rl  rl 

3T  J0uld^  +  2  J0QuuUld^  '  2  J0QuuXUl  d^  "-j06hd^  "  j06sd* 


(3.4) 


and  similarly  the  wave-energy  action  terms  from  (3.3)  become 


fr  L  uouidv  *  1 A  L  vi<5uUd  -  - 1  LviW  d* 


3  ^  3X  j  U1  u0  -  3a  gx  ^X^X  L^U0  ^  d'^ 

JU  V  v  u  v 


(3.5) 


=  "  j0ch  u0  d  v  '  I0‘3  U0d:;  » 

where  we  have  separated  the  terms  6.  resulting  from 

n  2 

the  variation  of  h  and  the  terms  resulting  from  the  0(c  )  perturbation 
and  the  slow  derivatives  of  Q,Uq,  and  k: 

5h  =  huUl  +  hu  *~XU1.  +  Vo^  +  hu  ^9XUlv^  +  20X<t>XUO'X^ 

A  A  A 

★ 

For  eh  to  be  a  dissipative  perturbation,  h(u,ux,uxx,X,T)  must  be  even  in  the 
u  argument  [i.e.  h(u,-u  ,u  ,x,T)  =  h(u,uv,uvv,X,T)].  Since  un  is  even  in  ip,  in 

A  XXX  XXX  U 

★ 

Generalizations  to  other  kinds  of  dissipation  (for  example,  dependent  on  third 
or  higher  derivatives)  are  easily  made. 


is  even. 


this  case  h(ug,exUg  ,exUg  ,X,T)  is  even  in  v.  hu  is  odd,  and  h^ 

-r  V  ^ 

We  thus  obtain  that  the  slowly  varying  terms  ha  e  no  effect, 


xx 


f  2 

f  d  .  =  1  un  d  .  =0,  if  ch  is  dissipative  and  e  h,  is  non-dissipati ve 

;0  s  jo  s  0  1 

(odd  in  the  ux  argument).  In  this  case  (3.2)  and  (3.3)  are  the  first  variation 

of  the  wave  action  equations.  Therefore,  the  wave  action  equations  (2.13)  and 

2 

(2.14)  are  valid  to  at  least  one  higher  order  if  ch  is  dissipative  and  e  h-j 

2 

is  non-dissipative.  This  means  that  (2.13)  and  (2.14)  are  valid  to  0(c  )  when 
u0  uo  +  eul  and  ~X  — <-  k  +  e  $x  »  if  the  contribution  from  6$  vanishes.  In 
general,  this  is  not  valid,  but  (3.2)  and  (3.3)  are  valid. 


4.  Modulated  Phase  Shift  for  KdV  Type  Equations. 

The  equations  for  the  modulations  of  the  phase  shift  q(X,T)  will  be 

determined  in  this  section  by  evaluating  the  appropriate  higher  order  action 

2 

equations,  (3.4)  and  (3.5).  We  will  assume  ch  is  dissipative  and  e  h,  non- 

r  1  ‘ 

dissipative  :  (l,un)d.  =  0.  These  equations  require  knowing  properties 
•  0 

of  u^,  which  satisfies  L(u-j)  =  R-j  •  By  differentiating  (2.3)  with  respect 
to  .,  we  immediately  obtain  the  well-known  result  that 


L(uq  )  =  0  ,  (4.1) 


that  is,  Ug  is  a  homogeneous  periodic  solution,  which  is  odd  in 

V 

It  is  less  well-known  that  the  other  homogeneous  solutions  may  be  obtained 
by  differentiation  with  respect  to  a  parameter.  Dobrokhotov  and  Maslov  [2] 
used  this  idea  for  both  ordinary  and  partial  differential  equations  of  the 
Klein-Gordon  type,  but  an  alternate  representation  introduced  by  Haberman  and 
Bourland  [3]  yields  a  better  result  for  c>(X,T),  independent  of  the  type  of 


15 


dissipation.  We  should  recall  that  the  parameter  ^  in  (2.3)  satisfies  the 
dispersion  relation  (2.6)  and  thus  is  a  function  of  A,  B,  and  k.  By  directly 
taking  the  first  partial  derivatives  of  (2.3)  with  respect  to  these,  we  obtain 


L(u0a)  =  ^Au0 

UuQg)  = 


(4.2) 


(4.3) 


Q,..,un  -  3a2k2uf 


JQ  '  "  "k  0  ’  ?wuu> 

k  V  'i> 


An  equivalent  form  of  the  latter  expression  is  obtained  using  (2.3) 


L(u0  )  =  (~k  "  k  ^  u0  +  QuuU0,  * 

k  y  'i> 


(4.4) 


This  has  yielded  a  second  homogeneous  solution,  ^gUQ  -  u^u^  ,  which  is  periodic 

A  B 

and  even  (and  hence  independent  of  u^  ).  To  obtain  a  third  homogeneous 
solution,  our  experience  with  the  Klein-Gordon  equation  suggests  we  should 
analyze  the  non-periodic  even  function  vyQ  •  By  direct  calculation 


L('Uq  )  '  ^uo  +  ?kQUuu0,  +  3i  k  U0 


Using  (2.3)  again,  we  obtain 


L(^un  )  =  2*un  -  kQ  u 


*  A  • 

uu  0 

V 


Thus  from  (4.2)  -  (4.4),  a  proper  combination  of  uQ  ,  Ug  ,  and  i|/u_  will  be 

^  A  k 

a  third  homogeneous  solution,  namely  wa(uq  +  £  'IWq  ^k  "  P  u0  *  which  15 

uk  v  UA 


•7i< 


16 


even  and  non-periodic.  [un  could  be  used  here  with  ..  -  „  instead  of  un  .] 

UB  B  °A 

We  are  now  able  to  solve  L ( u -j )  =  R^.  We  could  use  variation  of  parameters, 

However,  it  is  easier  for  us  to  note  that  from  (2.9) 


L(ui  >  =  Ri  ^  =  "Vo.  '  u  ■  Vo. 


L ( U 1  odd }  =  Rleven="U°l’  *  0<  V  Xx‘ "  (3Vo  x  +  3hXX"xV  }  +  h  • 


since  h  is  even  for  dissipative  perturbations .  The  integrals  in  (3.4)  and  (3.5' 


only  require  the  even  part  of  u ^ .  u -|  ,  significantly  simplifying  the  needed 

even 

calculations.  Equation  (2.3)  is  used  once  more  (to  eliminate  u^  )  so  that 


where  C^X.T)  is  an  arbitrary  multiple  of  the  periodic  and  even  homogeneous 
solution.  We  have  thus  shown  that  the  arbitrary  multiple  of  the  other  homogeneous 


solution,  Uq  ,  vanishes  when  the  secularity  conditions  equivalent  to  (3.4) 
and  (3.6)  are  applied.  Luke  [5]  had  shown  a  similar  result  for  the  Klein- 
Gordon  equation. 

The  substitution  of  (4.5a)  into  (3.4)  and  (3.6)  appears  formidable. 

We  can  make  the  calculation  elementary  by  using  an  observation  by  Haberman 
and  Bourland  [3]  for  the  Klein-Gordon  equation.  Equation  (4.5a)  is  equivalent 
to 


ui  =  3  •  (4.5b) 

ever,  J 

where  D  is  a  linear  first-order  (in  the  slow  parameters)  differential  operator 
with  variable  coefficients: 

0  :  :x  ^  Ta  +  C2(x’T)(-B  "X  '  “A  ‘  (4-5c) 

a 

In  s5,  we  show  that  D  is  the  Taylor  series  operator  in  the  parameters  A,  B,  and  k. 
This  operator  has  the  property  that 

A  A 

D ( k )  =  :x  and  D (,)=-: T  ,  (4.6) 

which  will  be  quite  useful. 

The  rhs  of  the  higher-order  action  equations  (3.4)  and  (3.5)  are  simplified 


by  noting  that 
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°  jnh^uo,vxuo,  *”XU0  ,X,T^  d  '  = 

*  ^  ♦'  ✓V 


f1  i 


J0  \  huDu0  +  \[~XDu0  +  D(kVJ  +  \/VV+  2^XD(k)u0 


'0 


'hd- 


and  similarly 


A  r1  2 

°  ,'o  h{u0*'xu0,  »‘x  U0  ,X’T)uOd'  =  *0'huOd'  ’ 


since  D  is  a  first-order  partial  differential  operator  in  the  slow  parameters 


that  commutes  with  the  fast  derivatives  We  also  used  u-|  =  D  uQ  and 

even 


D'k;  =  :  .  There  are  correspondi nq  simplifications  for  the  *hs  of  (3.4)  and 
'  3.  rj ) .  First  of  all, 


-1 


D  ,  Q..d 

'0 


-  fl 


rl 


U'  JnQUUDU0d’  ’  "  ;^QuXd-  =  ■  QmiXDund 


'0 


uuX  u  u0 


■1 


(unQ  -  Q)d  .  *  un  Q  Du.  d  .  ,  and 

,'q  u  u  ;0  0  Muu  0  ’ 


■1 


j  .0 


•0 


-  r,  id 
U  1  A 


0  uu  <  u  u0 


D  un  d  . 


since  0  depends  on  A,  B,  and  k  only  through  Up.  In  addition,  for  (3.5),  we  note 


a 


D 


d  . 


2k  D  (k)  ,  (u,,  )*■  d  .  , 

-0  d. 


where  again  D(k)  =  In  this  manner,  we  derive  the  remarkable  simple  result 
that  the  higher  order  wave  action  equations  (respectively  momentum  and  energy) 
have  the  following  nearly  conservation  form: 


*  'I  ,  A  '1  A  '1  }  -1 

-°;0U0d-  +?-°;0  Qud-  =  -  D;0hd-  -2^0QuXd- 


1  2  . 
2  u0  d 


2  'X 


D  ;()1(uoQu-0)d. 
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V^UA 
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1  *  1  3  2 

+  i  D  0  d  .  -  ■=  ; 

2  :0  X  2 


.  -  o  '1  o 

K  D(X  ;0  u0  d 


-  D 


uQd 


(4.8) 


in  the  case  where  the  perturbation  -h  is  dissipative  and  •  h-j  is  non-di  ss  i  pa  t  i  ve . 

It  is  interesting  tc  note  that  (4.7)  and  (4.8)  may  be  obtained  from  the  leading 
order  action  equations,  (2.13)  and  (2.14),  by  replacing  the  derivative  terms  of  t^e 
conservation  law  by  the  respective  derivative  acting  on  D  (i.e., 

as 

-y  by  -y  D  and  — y-  by  D)  and  the  non-derivative  terms  in  the  action  equations 
*  '*  -1  a  rl  -1  a  -1 

by  D  operating  on  them  (  i  ,e.  Q  d  .  by  0  ,  0„  d  ,  hd  ,  by  D  h  d  . ,  etc.  ) .  we 

•0  ;0  -0  ;0 

[3]  have  also  shown  this  phenomena  for  the  Klein-Gordon  equation,  and  thus  we 
believe  it  to  be  a  general  principle  for  appropriately  perturbed  systems.  Remap--- , 
a  general  proof  of  its  validity  could  be  given  using  Whitham's  [6]  averaged 
lagrangian  principle  with  the  necessary  inclusion  of  weak  dissipation. 


Although  (4.7)  and 


jrp  remarkably  simple,  expressing  the  validity  o‘ 
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-»st  .anation  of  each  wave  action: 


R  I  +  k  I 
DU  B  1T  k 


Al;<qA  +  BlxqB  +  klxqk  =  0 


(5.4) 


?  »e  note  that  ail  the  dissipation  terms  have  been  balanced,  leaving  a 

t  C-*  universal  form. 

*e  no*  e  the  conjectured  frequency  perturbation  (5.2)  to  eliminate  A-j 


(5.5) 


-  -  anq  k  =  :  .  It  can  be  shown  [non-tri vial ,  but  straight 

!  ’  d  ’  n a  4.5:  -  4.6)  ]  that 


,  ,  (5.6) 

C 


is  substituted  into  (5.4),  we  obtain 

T  .  :A  "  l-4C2^T  lB  +  ^BC2^TIA  +  *XT  !k 

H 

~  .  qA  "  (^AC2)XqB  +(wBC2^XqA  +  *XXqk  =  0  * 

M  X 

puj- valent  to  ■' 4 .  *  I '  -  (4.12).  In  this  way,  we  have  verified  the 
-.t-'-e  tnat  a  highly  accurate  dispersion  relation  may  be  obtained  using 
i.ec  : erturpat ions  of  the  amplitude  parameters.  Furthermore,  the  modulations 
•  '  •  .  ■ '  ■  *  t  are  dote'-  r*jd  from  equations  equivalent  to  the  highly 


accurate  wave  action  equations,  when  highly  accurate  dispersion  is  calculated. 

/V 

We  now  explain  our  earlier  remark  that  D  is  the  Taylor  series  operator  in 

the  parameters  A,  B,  and  k.  The  leading-order  oscillatory  wave  solution  for 

KdV  type  equations  is  in  the  form  Ug  =  Ug(yiA,B,k),  when  the  dispersion  relation 

(2.6)  is  satisfied.  Some  (but  not  all)  perturbations  of  this  can  be  obtained 

by  perturbing  the  parameters  A,  B,  and  k  (a  Taylor  series)  in  which  case 

u ^  =  D  uQ  ,  where 

even 


0  = 


(5.7) 


through  the  use  of  (5.5)  and  (5.6),  it  is  seen  that  (5.7)  is  identical  to  (4.5c) 


Appendix  A:  Secularity  Condition  at  0(e). 

2 

In  this  appendix,  we  show  that  eliminating  the  secular  terms  at  0(e  )  is 
equivalent  to  analyzing  the  0(e)  wave  action  equations.  From  (3.1a)  using 
(2.10),  we  immediately  obtain  (3.2)  using  the  periodicity  of  Ug.u^,  and  functions 
o"7  un.  Similar  results  are  obtained  for  wave-energy  action,  but  with  much  greater 


effort.  One  term  that  is  present  in 


R^Ugdv  =  0  can  be  evaluated  using  some 


ideas  that  worked  on  the  lOein-Gordon  equation  ([2], [3]): 
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2,2. 


C  L  ( u  -j )  “  “  &  ©xu  1  d  y;  ~ 

0  V  W'k 


J0 


u] R] . 


ave  integrated  by  parts  three  times  (equivalent  to  using  the  adjoint 

,  j'ed  the  definition  of  L,  and  noted  L(u-i)  =  R,,  where  R,  is  given 

1  rl 


J0 


R^Ugdi^  =  0  becomes  (3.3) 


af  ter  a  lengthy  but  elementary  calculation, 


Appendix  B:  Some  Integral  Identities 

The  averaged  densities  and  fluxes  in  (4.10)  and  (4.11)  can  be  related. 
By  integrating  (2.4)  using  the  periodicity  of  u  ,  we  obtain 


(B.l) 


a)  joV4<  +  2  k  j0Qud^  =  "  B  ‘ 


By  simply  using  (B.l),  (4.13)  is  obtained  from  (4.10).  It  is  not  much  harder 
to  derive  (4.14).  Multiplying  (2.4)  by  uQ  and  integrating  yields 

Ip  i  f  I  p  3  fl  p  r  1 

-  oj  u„  d  ^  +  j  k  unQ  dt^  -  a  k  un  di|is  -B  undi|/  .  (B.2) 

;0  u  c  ;0  u  u  J0  ip  J0  U 


However,  by  integrating  (2.5)  we  obtain 


u  2  d  p  +  k  [  Q  d  p  +  a2  k3  [  u  2  d  ip 

J0  'o  Up 


2  B  [  u_  d  +  2  A 
;0  u 


(B .  3) 


Equations  (B.2)  and  (B.3)  may  be  combined  to  yield 

3o2k2  jr1u2d;  =  ^  -  f  j1  u2  d  p  -  j1  (Q  -  uQ  )  d  p  , 

;0  Ip  K  K  Jo  Jo  u 

ri 

from  which  (4.11)  directly  becomes  (4.14)  since  terms  involving  (Q  -  uQ  )  d  ^ 

j0  u 

cancel . 
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